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CHAIN PROCESSES AND THEIR BIOPHYSICAL 
APPLICATIONS: 


PART I. GENERAL THEORY 


I. OPATOWSKI 
THE UNIVERSITY OF CHICAGO 


_ A mathematical theory applicable to the biological effects of radia- 
tions as chain processes is developed. The theory may be interpreted 
substantially as a “hit theory” involving the concepts of “sensitive vol- 
ume” or “target area”. The variability of the sensitivity of the organ- 
ism to the radiation and its capacity of reccvery between single hits is- 
taken into account. It is shown that in a continuous irradiation of a 
biological aggregate in which the effect of each single hit cannot be ob- 
served, recovery and variation of sensitivity are formally equivalent to 
each other so that a discrimination between these two phenomena is pos- 
sible only by discontinuous irradiation or by using different radiation 
intensities. Methods for the calculation of the “number of hits” and for 
the determination of the kinetics of the processes from “survival curves” 
or similar experimental data are given. The relation between the re- 
covery and the Bunsen-Roscoe law is discussed. The case in which the 
injury of the organism is dependent on the destruction of more than one 
“sensitive volume” is also considered. 


1. Introduction. The changes that living organisms undergo 
when subject to the action of a physical or chemical agent are often 
successive transformations: the organism, or a part of it, passes suc- 
cessively through various states before a certain condition accessible 
to a direct experimental observation is reached. A theory of such 
processes has a twofold objective: (1) to obtain some information 
about the states which cannot be observed in a direct manner, about 
their kinetics, their rates and their number; (2) to correlate the ob- 
served facts with the structure of the microorganism by determining 
that part of it which is the first biological cause in the chain of 
changes that the organism undergoes under the action of the external 
agent. That part, the sensitive volume of the organism (called also, 
perhaps in a less appropriate manner, its sensitive area), is formed, 
according to present viewpoints, of some heavy molecule or molecules 
which have a vital biological function and whose decomposition causes. 
some more or less permanent change of the organism, its injury, or 
death (cf. e.g. Jordan, 1939). Although obviously not all physico- 
chemical agents can be conveniently accounted for by such concept 
of sensitive volume, the latter seems to be suitable for an interpreta- 
tion of biological effects of radiations. 
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The usefulness of radiations as therapeutic means, their appro- 
priateness for investigations in various branches of biology, the very 
large experimental material which is available in these fields, and the 
hope that someday it may lead to a deeper understanding of the struc- 
ture of the living matter, create the need of a synthetic theory of 
these biophysical phenomena. 

The present work is developed from the existing mathematical 
approach in this field based substantially on the concept of a chain 
process. The paper gives practical means of application of the theory 
and shows-also some of its limitations. 


Since in a later part of the paper the theory will be applied to 
some experimental results on biological effects of radiations, it is con- 
venient to mention here the character and the order of magnitude of 
the phenomena in this field. The average minimum energy necessary 
to produce erythema of human skin by ultraviolet light is 1.8 X 10° 
erg/cm? at 2967 Angstrom (Blum, 1944, pp. 1145-1163). This corre- 
sponds to 10%" quanta per cell, if one assumes an average size of 
epithelium cell of 400? (cf. Hoerr, 1944, pp. 1385-139). Since the pas- 
sage of a photon into or through the cell can be identified, from a gen- 
eral viewpoint, with a new state of the cell, it is seen that we have to 
do here with an enormous number of states, which is further multi- 
plied by the fact that the capture of a photon by a cell is only a first 
step in a chain of transformations consisting of activation of mole- 
cules, ionizations, chemical reactions, ete. (cf. e.g. Blum, 1944, pp. 
1145-1163; Failla, 1936, pp. 87-122; Fano and Demerec, 1944, pp. 499- 
512; Lea, Haines, and Coulson, 1936; Noddack, 1933, p. 330). A simi- 
lar computation could be made also for irradiation by X rays, or to 
get a closer picture of the process, one couid identify the change of a 
state of the cell with a liberation of an electron in it or with a pene- 
tration of a primary or secondary electron into the cell. Theoretical 
formulae for such total number of electronic tracks in or through a 
sensitive volume of living matter were given by R. Glocker (1932) 
in relation to the wave length, the penetrating range of the electron 
and the size of the sensitive volume. In an irradiation of a sample of 
ascaris eggs by X rays, A. Zuppinger (1928, p. 755) estimated rough- 
ly an average number of 108 electrons liberated per cell at the time 
when one-half of the eggs of the sample were injured. These examples 
show how large is the number of events involved in an irradiation 
process. However, a mathematical analysis of experimental data re- 
veals sometimes a single event only; for instance, in many cases of 
killing bacteria by ultraviolet light or induction of gene mutations by 
X rays; in other cases the number of events is larger but seldom ex- 
ceeds an order of magnitude of 50 or 100. This apparent disagree- 
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ment is explained by the assumption that the sensitive volume is. ex- 
tremely small with respect to the size of the cell, so that the number 
of events accounting for the final effect is also small. 


2. General equations. To make the mathematics more concrete 
and in view of the applications that will be made in a later part of 
the paper, the theory will be developed by referring as an example to 
the biological effects of radiations. Consider a microorganism capable 
of assuming 7 + 1 states which will be numbered 0,1,---, 2. We 
will limit ourselves to the case in which an organism in a state i may 
pass directly only either to the state i + 1 (“direct process”) or to 
the state i — 1 (“reverse process”). In the irradiation of microor- 
ganisms the number of states is often interpreted as the number of 
collisions of photons with the sensitive volume, although there are 
possibilities of a much broader interpretation. In any case, the “di- 
rect process” is in this field very often a process of injury of the or- 
ganism and the “reverse process’, a process of recovery. The state 0 
represents then living organisms in normal conditions and the state 
m may represent dead organisms. Consider a homogeneous aggregate 
of N microorganisms and call NY;(t) the number of those micro- 
organisms which are in the state 7 at the time ¢. If we call ND; (t) 
the time rate with which the microorganisms pass from the state i—1 
into the state 7 (direct or injury process) and NR; (t) the time rate 
with which the microorganisms pass from the state 7+1 into the state 
i (reverse or recovery process), then 


dY;/dt=D;— Din + Ri —Rin. (1) 
The simplest intuitive assumption concerning the D;’s is 
D,=k Yiu, (2) 


where k; depends on the intensity of radiation and on the sensibility 
of the organism in the state i—1 to a given radiation. Consequently, 
the k;’s depend on the probability of collision between a photon and 
the sensitive volume. 

As to the rates of recovery, we will assume, following an idea of 
W. Swann (Swann and Del Rosario, 1931), 


Ri=9:1 Yin, (3) 


where g; represents the intensity of recovery of the microorganisms 
from the state i+1 to the state i, in a way similar as k; in equation 
(2) represents the intensity of injury of the microorganisms in the 
state i—1. The process of recovery does: notneed to be necessarily 
of a biological character because it occurs also in a variety of purely 
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physical phenomena; for instance, in the action of light on a photo- 
graphic emulsion it is known as “fading of activated emulsion grains” 
(cf. e.g. Silberstein, 1939; Long, Germann, and Blair, 1935). 

Consider the case in which at the beginning of irradiation 
(t = 0) all microorganisms are in the state 0. Then 


Y,(0) —1, Y;(0) =Ofori>0. (4) 


Very often we may look upon the state n as representing some effect 
of radiation; in many cases that state corresponds to the only observ- 
able degree of injury, although in some cases intermediate states (for 
instance, the degree of skin erythema) may be also discerned with 
sufficient accuracy for measurements or counting. The process of 
recovery, as represented by equation (3), does not involve the actual 
possibility of recovery from any state to all the preceding ones. If, 
for instance, a certain g;.. = 0, such possibility is precluded, and one 
could say that there is some kind of barrier between the states 7—1 
and 7 which the organisms cannot pass by recovering (see Sec. 6). 

Besides the idea of W. Swann, which is generalized in equation 
(8), other methods of taking into account the recovery are known in 
the literature. For instance, J. Reboul (1939) treats the recovery in 
a global manner as an exponential relationship between the number 
of totally recovered microorganisms and the time. Such approach 
does not take into account the different states or degrees of recovery. 
Other authors consider the radiation energy absorbed by the organ- 
isms and converted by them into a change of their state as decreas- 
ing exponentially with the time by a process of recovery (Williams, 
1934; Taylor and Pohle, 1938, p. 160; Holthusen and Braun, 1933, pp. 
204-205; Lea, 1938, pp. 493-494, 555-558; Lea, 1941, p. 609). 

Equation (1) involves the basic assumption that the process is 
sufficiently intense so that the representation of the time rate by the 
derivative dY;/dt is sufficiently accurate. Consequently, considering, 
for instance, the state n, the theory appears less accurate for small 
values of ¢, because of the conditions (4). 

The use of equation (3) involves the assumption that the time 
necessary for an organism to recover from a state to the preceding 
one is very short with respect to the total duration of radiation. A 
method taking into account the duration of recovery of each single 
organism will be discussed in a later part of the paper. 

We have the obvious relations: 


Do=ky=0; D;>0,h >0forlSizn, 
Dy 2 0, Ria 2 0, 
R=0, ghEOtorosisn=1, 61 hh 
Rn = gn=0 
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The last of these relations involves the assumption which we make 
throughout the paper that a recovery is impossible into the n-th state. 
Usually this is the final state of the process, and in such case we have 


Das = kn = Rn=Gn—=0, (6) 


The conditions kn. = gn = 0 may be interpreted also on the basis 
of equations (1) to (5), from which 


d »> Y;/dt ——— Kn Vi 1 In ae =: 0 
7=0 
may be easily derived, so that by equation (4), > Y; = 1, i.e., the to- 
i=0 
tal number of microorganisms in the states 0, 1,---, » is always the 
same at any time ¢. In other words, the states 0, 1, --- , n involve a 
closed system of transitions. 
If no recovery is possible from the final state, then also 


Rria=Gri=0. (7) 


For the sake of brevity, all equations will be written for a general 
state 7 and it is agreed throughout the paper that all quantities with 
negative subscripts are zero by definition. 

The intensities k; and g; will be assumed constant for each state 
2. Their possible dependence on time, as it occurs, for instance, dur- 
ing the mitotic division of a cell (cf. Henshaw, 1944; Pack and Quim- 
by, 1932) or in relation to its age (cf. Henshaw and Henshaw, 1933; 
Rentschier, Nagy, and Mouromseff, 1941; Sax and Swanson, 1941), 
may be taken into account by putting k; and g; in the form of a prod- 
duct of the type 


ks =¢(2) f(t), 9.=y(@) f(t), (8) 


and taking :) f (t) dt as the time variable instead of t (cf. Opatowski, 
1942). Equations (8) involve the assumption that either one of the 
two processes, injury or recovery, is not taking place, or that the ra- 
tio of the intensities of the two processes does not change with the 
time. ; 
; The fact that transitions between two consecutive states are con- 
sidered as the only possible in the direct process is a consequence of 
the assumption that the physico-chemical agent (e.g. the radiation 
quanta) as well as the biological aggregate are both homogeneous, 
that there is only one sensitive volume in each microorganism (cf. 
See. 7), and that the probability of a simultaneous action of two 
physico-chemical entities (e.g. quanta) in one sensitive volume is neg- 
ligible. The consideration of the same states during recovery may be 
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looked upon at least as an approximation. For instance, the photo- 
chemical processes of vision as they occur in the rods of human retina 
offer an example of such approximation; there are three states here 
which are determined by the products of a decomposition of the origi- 
nal components of the rods, the rhodopsin; the recovery occurs be- 
tween the second and the first state and directly between the third and 
the first, but the latter is of negligible amount (Moon and Spencer, 
1945). The possibility of obtaining biologically homogeneous aggre- 
gates through appropriate subcultures has been shown in the bacterio- 
logical field by H. C. Rentschler, et al (Rentschler, Nagy, and Mourom- 
seff, 1941; Rentschler and Nagy, 1942). 


3. Solution of the equations. A convenient method to determine 
the functions Y;(t) from equations (1) to (4) is to use the Laplace 
transformation ; that is, to associate with each function Y;(t) a fune- 
tion y;(s) defined by 


wie fees Y;(t)dt. (9) 


Then by known theorems (Doetsch, 1937, p. 153; Churchill, 1944, p. 
6), 


face (dY;/dt) dt = sy;(s) — Y;(0), 


so that multiplying both sides of equations (1) by e-*‘tdt and integrat- 
ing between 0 and co, we obtain the following linear system of equa- 
tions in y;’s 


(s + 2Io) Yo — Gos =1, 


| 
—kiyia + (8 + 21) ¥i— Gina =0 for 121 | ’ (10) 


where 27; = kis. + gi, so that I; is the arithmetic average of the in- 
tensities with which the microorganisms abandon the state i by direct 
and reverse process. 

The solutions of the system (10) are 


where 
Po=1,P,;=k,k-+-k; for 121, (12) 


D(n+1,n) =1, and any other D(i,n) is a determinant of order 
nm — i+ 1 whose only elements different from zero are those of the 
principal diagonal and of the two contiguous parallel lines, the He 
per” and the “lower” diagonals. 


I. OPATOWSKI 167 


Se ol pwns 0, 0, 
[Feiss , S te yy eee 5 Gir OF 
0 ke + 21, | 
Deny = ? b+2 9 s eed Geen Gir > . . . . . (18) 
Die ops Te Fe, Le ys fae kay or 
Ne Se eye vomras | ie Bos eet ae 


Determinants of this type are called continuants because of their con- 
nection with continuous fractions (Muir and Metzler, 1930, pp. 526, 
560, 561; Hellinger and Wall, 1943) ; they occur also in the theory of 
optical instruments (Rosin and Clark, 1941). An alternative form of 
D(i,n) which we necd is obtained from expression (13) by chang- 
ing the sign of each g and of each k in the upper and lower diago- 
nals. This does not change the value of the determinant by known 
theorems (Muir and Metzler, 1930, p. 25). The proof of equation 
(11) is obtained by applying Cramer’s rule. In fact, the determinant 
of the coefficients of the y;’s in the system (10) is D(0,n). The 
determinant A; obtained from D(0,%7) by changing its (i+1)-st col- 
umn into 1, 0--- 0 has the following character: if one draws a verti- 
cal line between the (i+1)-st and the (i+2)-nd columns and a hori- 
zontal line between the (7+1)-st and the (i+2)-nd rows, one obtains 
four rectangular arrays, such that the elements of the left bottom 
array are all zero. Such a determinant equals the product of the de- 
terminants formed by the left top and the right bottom arrays. The 
latter is D(i+1,n) and the former is easily seen to be P; (Muir and 
Metzler, 1930, pp. 103, 57), so that 
A;=P;D(i+ 1 xm); 

which proves equation (11). 

We need an expansion of D(i,n) as a polynomial of s. Call . 
D,(i,n) the value of D(i,n) for s =0, and M;(t,%”) a principal 
minor of order j of D, (i, n) ; that is, a minor whose principal diagonal 
is a part of the principal diagonal of D,(i, 7), then (Muir and Metz- 
ler, 1930, p. 72) 


Di +1) =s4[14+ 30% SoM +1,m) |, (14) 
j=1 , 


where Sy, stands for the sum of all the principal minors M; of 
D,(i + 1,7). In this way, each coefficient of s7 in equation (14) is 


‘a sum of (* =i , ) determinants of j-th order with the obvious relation 


M,i(i+1,n)=D,(i+1,%). 


For instance, for D(0, 7), i.e. fori + 1=0, we have 


168 CHAIN PROCESSES 


Tn 


Yu M,(0,n) =2351,, 
r=0 


r,S=n 21, 3 Or 
du M2(0,) = . (15) 
uM (0,7) 2 jen eed | (15) 
7,8,t=% | 7A by ’ gr ) | 
Due, (0,2) = 2 Kes 2hgeiee , ete. 
ee a POS Or | 


From expressions (11) and (14) the expansion of yi(S) in power 
series of 1/s may be calculated. In particular, for the n-th state we 
have 


Uns) = Ps/DO0..n) = Peis [1 ee One 


Su Mi(0,2)s? +---+D.(0,2)s**]} ? (16) 


y.(s) = Ps *"(1 + San, 5), (17) 


j=1 


where the a’s are obtained in terms of the Sw’s by the ordinary rules 
of division of power series, so that 


Ona =— D1 Ono = Dn? — D2} 
, (18) 
Ong = DP 2 Se Deke. | 
where for brevity 5; = Suv M;(0,”). By a known formula of the 
Laplace transform (Doetsch, 1937, p. 62; Widder, 1941. p. 94), we 
finally have from equation (17) 


Yn(2) =P,t| (n!)2 + San; t/(n+ 9)! ‘ (19) 
j=1 
which is obviously convergent for any finite ¢, because by equations 
(1), (2), and (3), Y,(¢) is a sum of exponential functions of ¢ multi- 
plied by polynomials of t. By equation (12), P, does not depend on 
_the g;’s, therefore in a sufficiently short experimentation (¢ small), 
the recovery process does not affect the function Y, (¢). 

Some conclusions of a general character may be derived from 
equation (19): In the field of radiations, the intensities of injury k; 
must be assumed, at least in a first approximation, as directly propor-: 
tional to the radiation energy per unit time, i.e. to radiation inten- 
sity J: 

k=Jki. (20) 
The reverse intensities g; , which represent the recovery process, may 
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be assumed in a first approximation as independent, or perhaps in- 
versely proportional to J since they may involve a physiological re- 
action of the whole organism. The well-known fact that tumor cells 
are not damaged by low radiation intensities (cf. Wood, 1944, pp. 88- 
89) but are destroyed at high intensities can be explained only by as- 
suming that higher intensities of radiation involve smaller intensities 
of recovery g;. The assumption of inverse proportionality between 
gi and J would involve an infinite intensity of recovery at zero inten- 
sity of radiation, which in itself cannot be considered as a contradic- 
tion of experimental facts. 

Now from equations (15) and (18) it is seen that a,,; is a homo- 
geneous function of degree 7 in the k;’s and g;’s. Therefore, if there 
is no recovery (g; = 0), the coefficient of ¢t’ in expansion (19) is of 
degree j also in J , by equation (20). Consequently, the effect of radi- 
ation, as represented by Y,(t), does not change if the intensity and 
the time of radiation change according to the condition Jt = const. 
We obtain in this way an interpretation from the viewpoint of the 
chain processes of the so-called reciprocity law asserting that a radia- 
tion effect depends on the total energy and not on its time distribu- 
tion (Bunsen and Roscoe, 1862). The reverse intensities g; take into 
account what is known in physico-chemistry as the failure of this law 
or the time effect; in biological fields they represent the recovery pro- 
cess on which the radiation-therapeutic concepts of fractional and cu- 
mulative dosis are based (cf. e.g. Taylor and Pohle, 1938, pp. 159-167; 
Holthusen and Braun, 1933, pp. 198-207). 


4. Numerical methods: First approximation methods. As has 
‘been mentioned, the most common case of application is the one in 
which of the n+1 states, 0 and m are the only observable states; ex- 
pansion (19) gives then a possibility of obtaining some information 
as to the laws governing the process from experimental data. In fact, 
from equation (19) we have 


log Y, —log P,,—logn!+nlogt + log| 1 + 1! On, U/(n + 9) 1] (21) 
Ea | 
For sufficiently small values of t, the sum > in [---] is negligible 
with respect to 1, so that the right hand side of equation (21) re- 
duces to its first three terms, and log Y, plotted against log ¢ is a 
straight line whose slope equals the number of states n and the inter- 
cept on the ordinate axis gives the geometric average of the inten- 
sities of injury P,1/". If larger values of ¢ are considered and conse- 
quently terms of second order in [ --- ] are taken into account, the sum 
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Ton 


Ana > Kren + Ora; (22) 
7r=0 
may also be obtained from experimental values of Y,(t). If the n-th 
state is the final one (kn. = 0), we may calculate approximately the 
arithmetic average G of the intensities of recovery Go, Gis °** Gna 
from equation (22), 


Ona ~ Pp" + NG, (23) 
where the arithmetic average of the k;’s is taken approximately equal 
to its geometric average. Formula (23) is exact if k, = k, =-+--= Kn. 


In all other cases it gives a value of G which is larger than the actual 
one because the geometric average is smaller than the arithmetic 
average (Hardy, Littlewood, and Polya, 1934, pp. 17-18). If the 
n-th state is the final one from which no recovery is possible 
(Kns1 = Yn = 0), then the above formula should be replaced by 


0.5 NPY ong G, (24) 


where G is now the arithmetic average of go, 91,°*- Yn2- 

The number of states n, which is sometimes interpreted as a. 
number of hits with photons, may be calculated also by means of the 
following formula, easily obtainable from equation (21): 


n= lim; [td log. Y,,(t) /dt] . (25) 


The use of formula (25) requires an accurate knowledge of Y, for 
smal! values of t, which experimentally is not easy to achieve for low 
radiation intensities J. However, since 7 by its physical meaning is 
independent of J, a sufficient degree of accuracy may be obtained if 
curves Y,(¢t) are known for various values of J. 


5. Formal equivalence between processes with and without recov-- 
ery. We will show in this section that.a chain process (p) with re- 
covery (g: # 0) between the states 0,1,---,m, of which 0 and n 
are the only observable states, is formally equivalent to a process: 
(p) without recovery (g; = 0) between the same number of states 
with appropriate intensities of injury kj. As to the process (p), we 
assume first that its n-th state ts final (kn, = 0), and that a recovery 
is impossible from that state (g,. = 0), so that J, = 0. We will see 
that the process (p) has also n as its final state, so that Kn..¥=0. 

Call Y,,(t; ki, gi) and Y,(t; ki, 0) the functions representing 
the n-th state of the processes (p) and (p) respectively. The proof 
of the above stated “equivalence theorem” consists in showing that: 
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for any set of values {k;, g;}, satisfying the conditions (5) to (7), 
there exists a set of n real positive numbers &; such that identically 
Y,(¢; ki, 9:) =Yn(t; i, 0). (26) 


Instead of proving the relation (26), it is more convenient to prove 
its Laplace transform, i.e., 


Males Ree Ui) bal Sts 5 0)s (27) 

or the following relations: . 
ky Se ee | are ee (28) 
D(0,n;ki,gi) =D(0,n;k;,0), (29) 


which by equations (11) and (12) imply equation (27). The symbol 
D(0, ) of equations (11) and (13) has been here replaced by the 
more explicit ones of condition (29), although for brevity D(0, 7) 
will be also used for D(0, 7; ki, gi). Since in the determinant ex- 
pression of D(0, n; &:, 0) in the form of (13) all elements above 
the main diagonal are zero, it is easy to see that, if the theorem is true, 


D(0,;k;,0) =s(s + ky) (st ke) --- (8 + Kn). 


Consequently, in order to show that there exists a set of n real posi- 


tive numbers k; satisfying the condition (29), it is sufficient to show 
that the roots of D(0,”,k;, gi)/s—0, ie. of (cf. relation (13) 


D(0,n—1)=0, (30) 


which is an equation of n-th degree in s [cf. relation (16)], are all 
real and negative. Then, these ~ roots taken with a changed sign are 
the values of the &;, --- , k, , because, as it will be shown, they fulfill 
also the condition (28). 

(a). The roots are all real, The determinant D(0,”—1) may 
be represented by the symbol ||a,,-||, where 


a,,—s+2I,, Giin=—9, Win=—ki 


for i= 0, 1,---, n —1, all other a’s being zero. The reality of the 
roots is proved by changing ||a,,.|| into a secular determinant ||D,,<|| . 
To do this, multiply each i-th row and divide each 7-th column of 
||a,,¢|| by a number x;, which will be now determined in a suitable 
manner. The new determinant ||b,,-||, whose value is obviously equal 
to ||a,,-||, is also a continuant with its main diagonal the same as that 
of the original determinant, 
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0:4 =a s ale 21; 5 
and the upper and lower diagonals respectively, 


Di iss = Fi in Ui/ Lier = Giki/Lin , (31) 
Di ia = Wiis Ui / Vi = Wii [Kia 


Now, if we choose x; in such a manner that 0ji4. = bisi1,i, then ||br,c|| 
is a secular determinant. Such a choice is possible by taking x; ac- 
cording to the relation easily obtainable from equations (31), 


Liv (9i/Kin)*/? Xi, 


which determines all x;’s by choosing e.g. 2. = 1. In this way, the 
secular form of ||b,,.|| and the reality of the roots in question is prov- 
en. It may be noted that in the new determinant 


Dijin = Dis, = (9i Kis) ee 


has a direct meaning because it represents the geometric average of 
the intensities of the processes occurring between the states i and 
basis 

(b). The roots are all negative. We prove this in the following 
manner: We show that all the principal minors of D,(0, » — 1) 
[ (cf. equations (15) and (16) ] are positive; consequently, D(0,n—1) 
considered as a polynomial of n-th degree in s [cf. the denominators 
in equation (16)] has all its coefficients positive. Therefore, by Des- 
cartes’ rule of signs, it has no positive or zero roots, and since all its 
roots are real, it must have all roots negative. It will appear from the 
following discussion that it is sufficient to consider three types of 
principal minors of D,.(0,n— 1): 

(i). Consider first the minor consisting of the first 7 rows and 
columns, i.e. 


ki, Jo> 0, 0, 

ki, ke + Qos 1, a. > - 
D,(0,7r—-1) = Oe; ke 5 ks + 91 Oongtss - ’ (32) 

Os ter ae eee ee a Pe 


and transform it by subtracting and by adding its rows in an alter- 
nating manner, i.e.: leave the first row as it is; subtract from the 
second row the first one; subtract from the third row the second row 
and add the first, ete. By known theorems we do not change the de- 
terminant by these transformations, but it is easy to see that the new 
determinant has all its elements below the principal diagonal equal 
to zero, and that the principal diagonal is k, , k, , ++. , k,.. Consequent- 
ly, 
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Dy(0.,7—1) Kika kk, > 0. (33) 


(ii). Consider now a principal minor obtained from D, (0;n—1) 
by omitting only a certain number of its first and of its last rows and 
columns. Such a principal minor is of the type [cf. equation (18) ] 


Kir + ir, Gi» 0, 0, 
: Kiss» Kise + gi, Jin» OS aria Sh 
Dy (% , m) =| 0 ? Kise ? Kiss of Jin» Ji+2 » ° ° , (34) 
0 7° . . . . ° . ° Kem ’ Kenaa ut Im-1 


where i 2 1, m S n—1. This determinant may be represented as a 
sum of two determinants having their first rows ki., 9:,0,0,--- and 
Jin,0,0,--- respectively, and all the remaining m— rows the same 
as in D.(i, m). The first of these determinants is of the type (32) 
and the second is easily seen to be gi.:D)(i+1, m), so that by equa- 
tion (383), 


Do(t,m) = Kiki +++ kma + gisDo(t+1,m). (35) 
We have directly from equation (34), 
D,(m, m) = kun ae Ina mL Be 


Consequently, the principle of mathematical induction applied to equa- 
tion (35) shows that D,(i,m) > 0. 

(iii). The third type of principal minor which must be aiacaasel 
is the one which is formed from D,(0, n—1) by omitting some of its 
intermediate rows and columns. Consider the case in which the 
omitted rows and columns are contiguous in the original determinant. 
Then it is easy to see that such a minor consists of four rectangular 
arrays, the elements of the left bottom array being all zero, and those 
of the left top and right bottom arrays forming each a determinant 
of the type (34) or (33). Therefore, the minor under consideration 
is a product of two determinants of the type (i) or (ii) and is con- 
sequently a positive quantity. 

It is clear from the discussion of the last type of minor that any 
other principal minor of D,(0, n—1) is a product of some principal 
minors of the above three types and is consequently a positive num- 
ber. In this way, the proof that the roots of equation (30), i.e. 
—k,,—Ke ,--- , —Kn are all real and negative is complete. But to ob- 
tain the equivatence theorem it is still necessary to show that between 
these roots and the original k;’s , the relation (28) holds. This is, how- 
ever, an immediate consequence of the fact that 


Ie, ley ++ Tig = [D (0 ,N—1) Je-0 = Do (0 nl), 
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which by equation (33) equals 4k, --- kn. 

The equivalence theorem which has just been proven shows that 
an analysis of experimental data on Y,(¢) cannot give any direct in- 
formation as to the recovery process if the intermediate states can- 
not be observed and no recovery occurs from the final state. Because 
of this latter condition the procedure of pure recovery consisting in 
obtaining experimental values of Y,(¢) after the irradiation has been 
stopped (cf. Fano and Marinelli, 1943) cannot be applied here either, 
and the only method to separate the intensities of injury k; from those 
of recovery g; is to carry out experiments with various intensities of 
radiation J and to take into account the different dependence of the 
k,’s and g;’s on J (cf. equation (20) and the corresponding paragraph 
of Sec. 3). 

As an example consider a process between three states 0, 1, 2, 
with 2 as the final state, from which no recovery occurs. By using the 
equivalence theorem, it is easy to see that such a process is equivalent, 
as far as the final state is concerned, to a process without recovery 
and with direct intensities: 


k, hes + (K*— kk) *? ) (36) 


where 2K = k, + k, + g.. The case k, = k, is the one considered by 
W. Swann and C. Del Rosario (1931) and it may be easily checked 
that their formulae can be obtained by using the relations (36) and 
the general equations of processes without recovery (Opatowski, 
1942). It may be seen from this example that the equivalence prop- 
erty as represented by relation (36) holds only for the resulting pro- 
cess between the initial and the final states but does not hold simul- 
taneously between intermediate states. The property of symmetry of 
the final effect Y,,(¢) with respect to the direct intensities of the com- 
ponent processes is also clear from relation (36). In fact, this rela- 
tion is symmetric in k, and k2 as well as in k, and ko . 

We have considered up to now in this section a process (p) in 
which knit = Gn1 = 0. If the n-th state is final and recovery occurs 
from that state, then ky. = 0 but gn. ~:0. A theorem similar to the 
equivalence theorem may be formulated also in this case, although it 
expresses only a proportionality of the effect. The theorem may be 
stated as follows: A number N of microorganisms which are initially 
in a state 0 and are subject to a chain process (p) between the states 
0,1,---,m involving direct and reverse transitions reaches the state 
m according to the same time relationship that an aggregate of a cer- 
tain number N of the same microorganisms would pass from the state 
0 to the state n if subject to a suitable chain process () involving 
direct transitions only between the same number of states. The differ- 
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ences in reasoning which are necessary in the present case to prove 
the theorem in question are the following: 

(a). Instead of D(0, 2; k:, gi) /s , we consider now the determi- 
nant D(0, »; ki, gi) with 27, = gn. [see equation (18)]. The ex- 
pansion of this determinant as a polynomial of s has all its coeffici- 
ents positive, except the term of degree 0 which equals D,(0,n) and 
consequently is zero by equation (33). Therefore, one of the roots of 
this polynomial is zero (Kn: = 0) and all the others —k, , —k, , --- —K, 
are negative, as before. The n-th state is final in the (#) process. 

(bo). The number of microorganisms which at t = 0 are in the 
state 0 of the process ()) must be taken equal to 


N=AP,/Pa)Ns.. where P,= tk, ky --- ky; 


and N is the initial number of microorganisms in the state 0 of the 
process (p). This change of the number of microorganisms is a con- 
sequence of the general expression of y,(s) [cf. equations (11), (12), 
and (28)]. To show that N # N, and consequently that we have in 
the present case only a proportionality of the effects and not an equi- 
valence of the processes (p) and (p), we have still to prove that 
P,,# P,. This can be done in the following manner: since one of the 
roots of D(0,n) is zero, P, , being the product of the non-zero roots 
of D(0, n), is the coefficient of s in the polynomial expansion of 
D(0, 7), i.e., by equations (14) and (16), 


P= du M,(0 ,). 


Now, all terms of this sum are non-negative. We show that one of 
them equals P, , and at least one of the remaining is positive, conse- 
quently P,, > P,. In fact, the term which is obtained from D,(0, 7) 
by omitting its last row and column is by equation (33) 


D,(0 , n—1) a 


The term which is obtained from D,(0, ) by omitting its n-th row 
and column is easily seen to be in the present case 


(Kn et In) D,(0 ’ n—2) —= k, ke ae Kaa » 0 . 


In this way the relation P,, > P,,, and consequently N < N, is proved. 
Since the process (#) occurs without recovery, it is intuitive that the 
initial number of its microorganisms in the state 0 must be smaller 
than in the process (p), which involves also recovery transitions. 

If the n-th state of the process (p) is not final and recovery oc- 
curs from that state, then kn #0, 9n1 #0. In this case all the n+1 
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roots —k,, —ky , --- —Kkn, of D(0,n) are negative. By equation (33) 
we have now 


D,(0,”) =f, bose = ee hae 


so that the above relation of proportionality of the effects can now be 
written in the simple form 


N= (Knia/ Kne N . 


In general, Kaa # kn, but if s = —Kyu is a root of the polynomial 
D(0,n; k;, 9:), then taking this root as —K,,, we have in the present 
case an actual equivalence theorem between the processes (p) and 
(p). This requires, however, a particular relation between the inten- 
sities k; and g; of the original process (p), which may be obtained by 
putting s = —k,,, in the equation D(0, ; ki, gi) =0. 

We can now discuss the question in what cases a complete transi- 
tion of the whole aggregate of microorganisms from the state 0 to the 
state n is possible. If the whole process occurs without recovery, 1.e. 
if all g;’s are zero, and if the n-th state is final (kn. = 0), we know 
that the transition tends to be complete as t > o because (cf. Opa- 
towski, 1942, p. 87) 


[Yn(t; ki, 0) Ji. =1. 


Therefore, a complete transition occurs also in the case in which re- 
covery is possible between the intermediate states of the process (p) 
and knw. = 9.1 —0. This is an immediate consequence of the equiva- 
lence theorem and of the equation (26). 

In the case in which ky. = 0 but gn. # 0, a complete transition 


does not occur because we now have for the number of organisms in 
the state n 


N [Ya (€: lees Gidlice = eee Olin = (PJE.JN.. 


It has been shown that P, < P,, therefore, not all microorganisms 
reach the state n even if t > o. This is intuitive because recovery 
occurs here also from the final state. 

If the complete process is without recovery and the n-th state is 
not final, the corresponding function Y,(t) tends to zero as ¢ tends 
to infinity (cf. e.g. Lundberg, 1940, pp. 49-64). This is true also if 
there is recovery and knu #0, gna ~ 0, because of the proportion- 
ality of the effects of processes (p) and (#). Therefore, the micro- 
organisms reach the state n up to a certain maximum number (which 
is necessarily smaller than N) and then decrease to zero as ¢ increases 
indefinitely. There is one maximum only and its location in processes 
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of this type has been studied by O. Lundberg and W. Feller (see Lund- 
berg, loc. cit.). 


6. Processes with recovery barriers. If recovery is impossible be- 


tween the states j and j — 1, that is, if g;. = 0 (cf. Sec. 2), then 
from equations (11) to (13) 


because in the determinant D(0, n) the elements of the first j YOWS 
and of the last n—j+1 columns are zero, so that D(0, ) equals the 
denominator in equation (37). The recovery barrier between the 
states j—1 and 7 splits the given process into two component pro- 
cesses: 

(1) the process which brings the organisms from the state 0 to 
the state j—1 with the direct and reverse intensities respectively equal 
tok, kz, wh kj; ; Jo, 915 °** Qj-2 and 

(2) the process which brings the organism from the state 7 to the 
state n with the intensities kj., , +++ kin3 Qj, *** Gna- 

Let NY;,,(t) be the number of organisms in the state n at the 
time ¢t, under the assumption that at the time t = 0, all the N or- 
ganisms were in the state 7. Let NY,j;.(t) have a similar meaning 
for the first component process. We have then from equations (37), 
(11), and (12) and by a known theorem of the Laplace transforma- 
tion (Doetsch, 1937, pp. 157, 163; Widder, 1941, pp. 91, 92; Churchill, 
1944, p. 37), 

Vom: == Vos ‘ aa , (38) 


where * is the known symbol of convolution defined by 
t 
g(t) * y(t) = exc wlt—1) dr. (39) 


In this section Y,,; is used for clarity instead of the previously used 
symbol Y;. 

All formulae of the previous sections are still valid in the pres- 
ent case, but if G means now the arithmetic average of those inten- 
sities of recovery which are presumed to be different from zero, then 
the term nG should be changed in equations (23) and (24) into 

n—1)G. 

: tee discussion of this section may be easily generalized to the 
case when recovery barriers exist between several states. One ob- 
tains instead of equation (38) the Ree formula: 


Yin a Yoja * bere m-1 eS "ga 1 Von, (40) 
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which appears as a generalization of a relation valid in the case in 
which the whole process goes on without recovery, that is, when re- 
covery barriers exist between each two states (Opatowski, 1942, p. 
85). 


7. Organisms with more than one sensitive volume. To treat this 
case it is convenient to interpret the function Y;i(t) = Yoi(t) as a 
probability that a microorganism be at the time ¢ in the state 7 if it 
is at the time t = 0 in the state 0 (cf. Opatowski, 1942, p. 84). Con- 
sider for simplicity the case of two sensitive volumes (v) and (w), 
and discuss first the case in which a certain effect (€) is reached in 
the organism when (v) is brought from a state 0 to a state p, and 
(w) from another state 0 to a state g. Assume that the changing of 
a state of either one of the two sensitive volumes is independent of 
the state of the other sensitive volume. Since a state i of the micro- 
organism may be characterized as a simultaneous presence of (v) in 
a certain state p’ and of (w) in a certain state q’, it is clear that a 
microorganism may reach the effect (e) through different intermedi- 
ate states (cf. Sec. 2), so that the theory developed in the previous 
sections cannot be applied without modification. This, however, can 
be easily obtained by probabilistic considerations. In fact, the prob- 
ability Yo.,,~(¢) that the microorganism shows at the time t the effect 
(e) if at the time t = 0 it was in the state 0: that is, if (v) and (w) 
were both in the state 0, is 


Yo;0,q—= Yor Yog- (41) 


Equation (41) is substantially a formula of H. Bethe and W. Seyfarth 
(cf. Glocker, 1932). Instead of relation (25), we have in the pres- 
ent case 


p+ q=limeso [td loge Yo.,¢(t) / dt] . (42) 


An application of the power series expansion (21) gives now 


log Yo;»,q = log Pp. —log (p!q!) + (p + q) logt 


+ log (1 + at +...) , eae! 


where a = (p + 1)"a,,. + (q +1)7a,,, and a,,,, and a,,, have a simi- 
lar meaning as a,,, defined by equation (22). 

The case in which the effect (e) is reached when (v) reaches the 
state p or (w) reaches the state q may be treated in a similar manner 
on the basis of general theorems of probability. 
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THE FLOW OF A VISCOUS FLUID IN AN ELASTIC TUBE: 
A MODEL OF THE FEMORAL ARTERY 
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The solutions of the classical differential equations for flow of a vis- 
cous fluid in an elastic tube are shown to yield a description of some 
aspects of the flow of blood in the femoral artery. The parameters in- 
fluencing the speed of the pulse wave, the pressure, and the velocity for 
non-pulsating and pulsating flow are exhibited. The relation between 
pressure and velocity are considered in some detail. It is shown that for 
pulsating flow the slope of the pressure versus fluid velocity curve is 
equal to the pulse wave speed multiplied by the density of the blood. 


The mathematical treatment of blood circulation has the usual 
intractable features of general hydrodynamics as mentioned by N. 
Rashevsky (1945a, b) further complicated by such factors as the 
variation of the viscosity and of the compressibility of the blood main- 
ly due to the corpuscular structure, as well as by the nervous and 
the chemical control of the size of the blood vessels. In spite of the 
present impossibility of even a good approximate general treatment, 
the classical equations do yield results which have significance for the 
physiologist and biophysicist. 

This paper discusses the equation of flow of a viscous fluid in an 
elastic tube and uses data on flow in the femoral artery to test the 
solutions. 


The Differential Equations of Flow: 

The equations for this problem may be derived by the usual meth- 
ods as exemplified by O. Frank (1926) or A. Foch (1932). The equa- 
tions are 


ou Oop 
—+—=yV" 4, (1) 
Be ERE ee 
ou 10 
apatite A (2) 
Cre. voat 


where 


/ 
w= velocity along the axis of the tube 
(In x-direction) 
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‘= pressure 

p = density of fluid 

w= viscosity of fluid 

E=} the modulus of elasticity of the tube wall 
VY2u= the Laplacian of u 

K= the compressibility of the blood 

e= thickness of the tube wall 

To= radius of the unstretched tube. 


In our treatment, p is constant and, therefore, K = 0. 
Differentiating equation (1) with respect to t and equation (2) with 
respect to x, and eliminating the p-term, we find 


On expressing the Laplacian in cylindrical coordinates and substitut- 
ing u—=C exp (—at— 62x) R(r), we are led to Bessel’s equation for 
R(r) (Von Karman and Biot, 1940), 


@k- 1 -dk 


——+/7R=0, 
d r? radar 
where 
pas?—po? + (na—pv?) B. (3) 
The final solution is then 
u=exp(—at—fx)[C,J.(Ar) +C.Y.r)]. (4) 


Inasmuch as u must be finite at the axis and since 
ey lie) > —o as) oe 0 OS 


the complete solution is therefore 


U=ZTCuIo(Anr) exp(— an t — Bn x). (5) 
In our flow problem, we are primarily interested in the pulse 
wave speed, the pressure, and the fluid velocity. 


The speed of the pulse wave: 

If we are treating flow of a non-viscous fluid, the u-term on the 
right of equation (1) would be absent. Elimination of the u-term be- 
tween equations (1) and (2) would lead to the wave equation for p 
where the speed of propagation of the pressure wave, v, would be 
given by 


Bas 
v= FI 
\ 2 po 


(6) 
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This is the well-known Moens’ equation (Moens, 1878). Although de- 
rived for a non-viscous liquid, it has been found to apply quite satis- 
factorily to viscous flow. 
In the v-equation we seem to have a valuable equation for finding 
a relation among some physically significant parameters (e, E , p.»To)- 
These parameters are, however, most difficult to measure and even to 
define satisfactorily. J. C. Bramwell and A. V. Hill (1922) found that 
Moen’s equation held for flow in an artery but they transformed it by 
replacing the original parameters by the compressibility of the artery 
which is more easily measured: 
v = 0.357 J V Bele ; 
adv 


> 


where p is measured in mm of mercury (Hg), V is the volume per 
unit length of the artery, and v is in metets/sec. 


Non-pulsating flow: 

The elastic tube is connected to a reservoir where the pressure 
may be maintained constant for any desired length of time (i.e., at 
x=—0,p=—p,, for all t’s). The tube will stretch to some maximum 
radius, b. At the tube wall (r = b) we assume the usual Poiseuille 
condition, u — 0, hence 


Jo(Ab) =0. ’ (7) 
Using equation (3) and equations (2), we get for the pressure 
v? pp 
p= (x)— rere (8) 
At the wallr =b,u=—0O and 
f(x) =D. 


But also at wall, the pressure is distributed in overcoming the ten- 
sion of the walls and the outside pressure, so 


p=p,ten(>——), (9) 


where p, is the pressure in the medium surrounding the tube. If b 
varies along the tube, as it probably will, », is not constant. For vari- 
able b there would have to be a different value of 4 for each point. 
A more tractable procedure would be to treat the tube as composed 
of sections of constant radii, but the radii decreasing as we go from 
section to section. For each section then there would be a series of 
solutions corresponding to equation (7) for the radius of that sec- 
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tion. There is, however, no need for our doing that for this discussion. 

The data of S. R. F. Whittaker and F. R. Winton (19383) may be 
used to test equations (3) and (8). Whittaker and Winton (1933) 
perfused the hind leg of a dog with whole blood through the femoral 
artery and recorded the relation between pressure and velocity of 
flow. The flow took place through the entire complex of artery, con- 
necting arterioles, capillaries, and finally into the veins which were 
cut to allow the liquid to run into a collector. The equations used in 
our derivation are for a single vessel; hence, in all correctness, we 
may consider our theory as yielding the single vessel equivalent of 
their complex system. 

Their data may be expressed approximately as 


p=20+44, (10) 


where j, the average pressure in the cross-section of the artery, is 
in mm of H, and @ is in cc/sec while our p is in dynes/cm? and u is 
in cm/sec. Using equation (4), we have for the flow in cc/sec, 


6 (0) b 
u=2n i! rudr = 5 2.2C, exp(—a,t— a2) f Tid Ager) adr ; 


0 n=1 0 
since 


d 
— [a Jae) a= ee tea (SE). 
dx 

We have corresponding to equation (5) 


“=22b > Cy dnt Ja (dnb) exp(— a, t— B, 2) , (11) 

where the J,,(2) is the Bessel function of the first kind of order m 
(Von Karman and Biot, 1940). 
Equation (11) should be able to satisfy the conditions here imposed 
and the pressure equation derived from it should satisfy at least ap-. 
proximately the condition p = p, at x = 0 for allt. 

The relation between and @ is of the same type as that between 
p and wu, that is, equation (8). We see then that a and 8 may both 
be real, that 6 is negative and that 


a v? p B a pa ) 
a 2\ nbz /* 
The preceding equation relates a and 6. From the roots of equation 


Al ya A = 2.40/b; A, = 5.52/b, --- , -»» we have the various J’s to be 
substituted into equation (3) so that we may determine both a and’ 
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B in terms of p,v,u,and b. The algebraic expressions are easily ob- 
tained and show that there is a choice of sign so that a > 0 and that 
a varies inversely as yu. 

Although the preceding analysis seems quite adequate for a de- 
Scription of the flow in the hind leg of a dog, so far as the relation be- 
tween pressure and velocity for non-pulsating flow, the data of T. E. 
Machella (1933) reveal that p and w are also linearly related for the 
femoral artery proper. His data were for pulsating flow, which we 
shall consider in the next section. 


Pulsating Flow: 

In pulsating flow, the pressure and velocity are some periodic 
functions of time and possibly also of position. We assume that ve- 
locity, wu, and pressure, p , may be represented as a Fourier series of 
time. In order that our solution may be amenable to such expression, 


we may take 
a=a+t@m and p=l+im. 


Equation (3) becomes on equating real and imaginary parts 
ps? @= p(a? — w*) + (ua— pv?) (P—m) —2ulmk, 


12 
ui wo=2paat+2lm(ua—pv?) +uo(?—m’). (38) 


The velocity must be zero at the wall as before; thus 4 will again satis- 
fy equation (7). We still have four parameters (a, w,/, m) and 
two equations (12) besides the boundary conditions. 

T. E. Machella (1936) gave curves showing the relation between 
the velocity and pressure curves at the same point in the femoral 
artery (see his Figure 7). The curves have the same shape, but the 
numerical values of maxima and minima differ. We shall show that 


solution for u fits his data. We may write 


i > F(a, r) exp(—a, t) [C, cos a, t + C,’ sin wy f] . 
n=1 
Since x is fixed and r may be considered as constant, we must ascer- 
tain whether a u of the form 


u=> exp (— Oy t) [C,, cos w, = + C,, sin w, t] 

n=1 ; : 
can fit Machella’s results for both u and p. We have a choice. Hither 
we may fit our uw for a single cycle of systole and diastole or fit a 
complete curve of many similar cycles. If we choose to find u for the 
general curve with many cycles, a, —0, and our problem becomes one 
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of finding the C, and C,' to fit the experimental data by harmonic 
analysis. The author reproduced Machella’s curve for u and p found 
the coefficients by a numerical method (Von Karman and Biot, 1940). 
By taking a sufficiently large number of terms the series may be made 
to fit as closely as desired. 
The expression for p is from equation (8) 
as2 * 
eS wa (13) 


@ 


Thus if the velocity, wu, can be expressed as a Fourier series with real 
coefficients, the pressure, p, must be expressed by the same series. 
with each coefficient multiplied by a constant factor. Machella’s p 
ranges from 88 to 156 mm of Hg and his wu from 5 to 93 cm/sec. These: 
two points are sufficient, of course, to determine the slope and inter- 
cept in equation (18). However, equations (1) and (2) are not ap- 
plicable for the assumption that a, = 0 involves a non-viscous fluid,. 
for in order to satisfy equation (3) and the experimental conditions, 
it is necessary that u = 1= 0 which is the non-viscous flow problem. 
The relation between and w is still given by equation (13). But the 
differential equation now requires that 


o=vmM, 
so that equation (10) becomes 
p= (vp)ut pM. (14) 
Machella’s (1936) curve 7 gives 
p = 0.773 u + 84.18, 


where # is in mm of Hg. In order to compare with equation (14) we 
must multiply by (1,018,000/760) to convert » to dynes/em?. We 
have then 


. p=1020u + 110,000. 


_ Therefore, the Fourier Series calculated for u gives p through this 
linear equation. Probably the most interesting aspect is that the slope 
of the p—u curve gives the speed of the pulse wave, for since p © 1, 
the slope is very nearly v. In Machella’s (1936) experiment we have 
that v = 10 meters/sec which is a reasonable value. Alternately know- 
ing v, we can calculate the compressibility or the other pete sit 
from. Moens’ equation. 

The preceding analysis is not quite as attractive as it seems, fer 
the pulse velocity has been found to vary with the pressure. J. C. 
Bramwell and A. V. Hill’s (1922) values for pulse velocity as a func- 
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tion of effective pressure, (p — po), of an excised section of artery 
reveals that as (p — »)) varies from 0 to 70 mm of Hg, the pulse 
velocity goes from 1.5 to 4.9 meters/sec. This experiment was for 
non-pulsating flow and is relevant to our equation (9). We cannot 
explain the variation, however, by a mere increase in radius, since 
the velocity for the unstretched tube, is too small (v, = 1.5 meters/ 
Sec). The logical explanation is that the elasticity, #, varies with 
pressure. Equation (9) is still valid physically; and it can be used 
to determine the variation of E with effective pressure by measuring 
the radius. Equation (14) is still valuable in that it yields an ap- 
proximate value of the pulse velocity. 

If we are interested in finding w for one cycle only, we may take 
a, ~ 0, and we have a slightly more difficult problem in harmonic 
analysis. From Machella’s (1936) data, we may determine the co- 
_efficients of wu. In place of equation (13) we have 


v? p(al+mo) 
See ee 


‘es C= ae ae 


In order that p may be real the imaginary part of » must vanish: 
lwo=—am. (15) 


The coefficient of u here must equal 1020 gm/sec cm? as before. The 
equation resulting from equating the coefficient to 1020 gm/sec cm? 
with the equations (5), (7) and (15) determine the relations among 
the parameters, a,1,m,A (since w is determined by the period of the 
cycle), and p, v, and u. 

Although the solution for a single cycle may initially seem rather 
artificial, the fact that it involves the complete equation recommends 
it to the author. We look upon the flow as being initiated at each 
systole, continuing the diastole and the values of u and » being altered 
by the characteristics of the elastic tube (the femoral artery). 

It is interesting to note that Machella’s pressure—velocity curves 
for the carotid artery do not show similarity of shape. This lack of 
similarity may be explained by the nearness of the carotid to the 
heart, thus the surge is so great that the approximations inherent in 
our differential equations are much too crude. 

The author wishes to thank Professor N. Rashevsky for the op- 
portunity of examining his recent papers on flow (1945a, b) in manu- 
script. He is indebted to certain members of the Howard University 
Medical School for data. 

The publication of this paper was assisted by the publication fund 
of the Natural Sciences in the Graduate School of Howard University. 
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; Under the assumption that the elongated shape assumed by a grow- 
ing nerve cell is caused by surface tension forces between the cell, its 
fluid medium, and a fibrous substrate track along which the cell grows, 
equations of elongation and conditions insuring elongation are derived. 
One specific type of cell-substrate contact is considered. Equations of 
elongation and conditions for elongation are treated in both the non- 
frictional and frictional types of motion of the cell-process. 


In this paper we propose to develop a theory of the growth of 
nerve fibers. We shall explore the consequences of various assumed 
relations between the physical factors such as interfacial tensions, 
cell volumes and shapes, etc., which might plausibly enter as deter- 
mining this growth. 

At this point of theoretical investigation, we are chiefly inter- 
ested in the role played by these physical factors when related in cer- 
tain ways. In order to discover the influence any one of these factors 
may exert, a thorough study of many cases must be attempted. 

Certain conclusions will be forthcoming which may aid in en- 
abling us to decide the basic question of the mechanism of the growth 
of nerve cells, both in the living organism and in the tissue culture. 

The growth of nerve fibers in tissue cultures has been described 
by P. Weiss (1934, 1941) as a drawing out of elongated threads, called 
filopods, from the central body of the cell. 

These elongated threads grow out along fiber tracks which are 
produced from an exudate originating from the cells in the culture. 

We shall discuss the process of growth of a filopod under the fol- 
lowing conditions: (1) the growing cell-process may be considered as 
a viscous thread approximately cylindrical in shape in contact with 
the substrate fiber; (2) interfacial tension forces at the interfaces 
of cell-process, culture medium, and substrate fiber are predominant 
in causing elongation of the cell-process; (3) the volume of the cell- 
process remains constant during elongation. 
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I 
We shall first suppose that the cell-process makes contact with 
the substrate fiber along a thin strip. Figure 1 shows both a longi- 


eee ee 5 ee ne 


Cell 


Cell Process 
Process 


tudinal and cross-section of the cell-process and substrate in this case. 
We assume that the radius of the substrate fiber is much greater than 
the radius of the cell-process. Hence we may consider the surface of 
contact between cell-process and substrate approximately as a fiat 
rectangle. 

We denote by” 


= width of the cell-substrate interface 

7 = radius of cell-process 

z= length of cell-process 

6 = angle subtended by p 

B= angle of contact of cell-process with substrate. 


It has been suggested that the filopod grows because it is “pulled 
in” between the substrate and the medium, in a process similar to 
wettening. But in such a case an isotropic interfacial tension would 
result in the spreading of the cell material in a thin layer all over the 
interface between medium and substrate. A growth in one direction 
is possible only if the interfacial tension at least between two of the 
three component materials is anisotropic. Inasmuch as the substrate 
is supposed to consist of oriented long-chain molecules, this unavoid- 
able anistropy is made physically plausible. In this case, the inter- 
facial tension becomes a two-dimensional tensor. Approximately we 
shall treat the problem as if we had simply two distinct interfacial 
tensions in two perpendicular directions. 


Let us denote by I’,, the interfacial tension of the cell-substrate 
(cs) interface parallel to the z-direction and by y.. the interfacial 
tension of the cs interface normal to the z-direction. 

We also denote by I'm. the interfacial tension of the medium-sub- 
strate (ms) interface parallel to the z-direction and by yom the inter- 
facial tension of the cell medium (cm) interface in any direction. 


RICHARD RUNGE 191 


Let us denote by 


S-s = area of cs interface = zp 
Sms = area of ms interface ee (1) 
Sem = area of cm interface = 4r2(2n — 0) + r(2n—6)z 


From elementary geometry we have 6 = 28, and since f is an angle 
of contact, we may assume that it, and hence 6, are constant. Thus 
from expression (1) we have 


dS cg = zdp + pdz, (2) 
and 
dS em = 2(a— B) [rdz + (z+ 7r)dr]. (3) 


In equation (3), z >> rand dz >> dr, hence we may ignore the term 
rdr and obtain 


dS om = 2 (a — B) (rdz + 2dr). (4) 
We also have 
GSms = — AS., = — zdp — pdz. (5) 
From Figure 1 we see that 
p—2rsin B, (6) 
thus equation (2) becomes 
dScs = 2r sin B dz + 2z sin fdr. (7) 


In equation (7) the part 2r sin § dz represents an increment of S,s in 
the z-direction while 2z sin 8 dr represents an increment of S., nor- 
mal to z. 

If we denote by dE the change in energy resulting from a change 
in length dz , we have from equations (4), (5), (6) and (7): 


dE = 2r[sin 6 (Ics — F'me) a8 Vou (te P) 1 dz + 


2z[sin B( yes —I'ms) + ¥em(a— B)] dr. Sh 
The relation between z and 7 is given by 
7 (a — B)z = V = volume of cell process, (9) 
and since V is assumed constant, we have 
spe lh alaig (10) 
eel 2) Ee 


Combining equations (8), (9), and (10) gives 
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; Vdr c 
= i — L'ms) 'F Yom ok Fact eg a 
dE A4[sin B (Ics — Ems) ‘+ Yom (a — B) Ge 
(11) 
oY — [ain B (ye Dna) + You = 8) oF 
Te | bd! cs” £ms em \Jt — . 
1? (a — B) id Y 
If the process is to elongate spontaneously, its elongation must 
be accompanied by a decrease in energy. Hence 
dE <0. (12) 


Since V > 0,7? > 0, and # < a, equation (11) implies that inequality 
(12) is equivalent to 


{—2fsin f ( 6, Ime) + (xa — B) coon | +r 
[sin B (yes — on) ok (2 — B) veal} dr <0. 


In order to have elongation, we must have dz > 0, or, because 
of equation (10), dr < 0. Thus in order that inequality (12) holds, 
we must have: 


sin B (ys m9 Lee ee) - (= P) Fou a O.. (14) 


(13) 


Since inequality (14) is independent of r and z, we may conclude 
that once elongation begins it will continue until the end of the sub- 
strate is reached. 

This conclusion is limited by the fact that as 7 approaches a mag- 
nitude comparable to molecular dimensions the preceding analysis of 
surface tensions is inadequate and equation (8) does not hold. 

We shall now derive an equation giving z as a function of the 
time t. We must at first content ourselves with such an equation of 
elongation as can be derived assuming that z = z, when t = 0, where 


2 itself represents an elongation sufficient to make the cell-process 
cylindrical. 


As a first approximation, we may consider the filopod as a vis- 
cous thread, pulled by a force concentrated at its end and equal to 
—dE/dz, choosing the direction of elongation as positive. Denoting 
by 7 the viscosity of the thread and by A, and z its cross-sectional 
area and length respectively, we have (Nadai, 1928) 


1dz 1 —dE/dz 
PW ole PAS 
In our case, Ay is given by 


(15) 


A, =r (a— B), (16) 
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while —dE/dz is given by equation (8) as 


—dE 
a, orlsin B(L'ce — I'ma) + Yom (a — 8)] — 
(17) 
dr 
22[sin B(ycs — F'ms) + Yom (a — B)] —. 
dz 
Combining equations (17), (10), and (9), we obtain 
dE ; 
Fe lve + Dns — 2P'en) sin B— (A—B) yom]. (18) 


Combining equations (18), (16), (15),.(10), and (9), we have 


dr 
—67 (x — 8) gg SE eet Lm — 20 ce) —=" (ct 5) von (19) 
Putting 
sin B (yes 2 Poe 20) SF (x — 8) yom 
——————— Fhe 7); (20) 
—6n (x — 8) 
we have 
LS (21) 
din 
Integrating equation (21) gives 
r=Qt+7r. (22) 


From inequality (14) and equation (20), we see that Q. < 0. From 
equation (22) it follows that after t = —7,/Q, seconds r = 0 (and 
hence, z= ). Actually, however, as r approaches molecular dimen- 
sions, equation (22) does not hold. 
From equations (9) and (22) we obtain 
= Sees Aer raey (23) 
(a— 8) (To + Qot)? 


In Figure 2 a graph of z is drawn according to equation (23). The 
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dotted line represents equation (23), whereas the solid line takes into 
account the fact that elongation ceases as r approaches molecular 
dimensions. 

In the preceding discussion, the possible effect of all frictional 
forces was not considered. The internal friction present in an elon- 
gating viscous body is taken into account by equation (15). 

The coefficients I’., and y.; measure the net amount of attractive 
force between the molecules of substrate and cell. As the cell-process 
elongates, the number of molecules within the cell which are brought 
sufficiently near the substrate for molecular attractive forces to affect 
them increases. Thus the total net attractive force between cell- 
process and substrate increases proportionally to Ses . 

We may take into account this change by considering the net 
attractive force between substrate and cell. as a normal force acting 
to produce a frictional drag force upon the cell-process. Since the co- 
efficient of static friction u., is given by (Foley, 1941) 


Mes = f/F n, (24) 


where F', is the net attractive force between cell and substrate, and 
where f is the frictional force, we may write 


f = LcsKesS es : (25) 


In equation (25) k.. is a constant of proportionality. Assuming that 
the coefficient of static friction is constant and putting 


fo aS PcalGes , (26) 
then from equations (25), (1), and (6), we have 
f=2f er sin B. (27) 


The force F" stretching the thread is then given by equations (27) 
and (18) as 


F’=r7[ (yes ~f | get —~ BF'ca) sin p-— (1 — B) yom — 2fo zsin 6). (28) 


The expression (28) now takes the place of —dE/dz in equation (15). 
Combining equations (15), (16), and (28) and putting 


(vce FD 2 i a 4 sy 
Ric Bah Aiand Bacal ail A WEG |S 
3nV (xn— B)V 
and 
i =), (30) 


8nV(a—B)V 
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we obtain 
dz 
as tt zy (31) 
Equation (31) is in a convenient form for purposes of comparing the- 
ory with experiment. Setting dz/dt = 0 in equation (31), we find as 
a final value for z 
ee Ry il (Ves + oe a Pd Bre | sin f a (2 — B) Yom 
oo - 
So 2f, sin p : 


Since elongation, even when frictional drag is not considered, is pos- 
sible only if equation (14) holds, we may state that equation (14) 
must hold also in the present case for elongation to occur, and there- 
fore, — R,./S, > 0. Hence if 
tore: 0 
ea p. 2 
3 (33) 


0 


(32) 


then dz/dt > 0 by equation (31) and z will increase until z= — R,/S,. 
The same result for a final value of z is obtained by making use of 
equation (28) and 


ey (34) 


as well as the conditions dE < 0,dz>0. 
Integration of equation (31) gives 


2 V—-S.. Rot+tvVzvV—RSo | 


RoVz R.V Ro R, al Vz V—E.S> 


where 
2 =, Ro + V—ZcltoSo mt 
uh SS pe (36) 
Rov 20 RVR, Ten Pr V—Zl So 


h,=— 


At 
A different type of cell-substrate interface arises when the cell- 
process, acted upon by isotropic interfacial tensions, surrounds the 
substrate fiber completely. Denoting by a the radius of the substrate 
fiber and by r the radius, measured from the center of the substrate 
fiber, of the cell-process which may be approximated by a cylindrical 
sheath covering the substrate, we have 


Sem = 2arz'+ n(r?— a?), (37) 
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and 
Ses = 2naz . (38) 
Hence 
dS om = 2ardz + 2nzdr + 2ardr , (39) 
and 
; GScs = — TSne — 2nldZ . (40) 


In equation (39) we may drop the insignificant term 2z7dr since 
z>>~r, obtaining 
OScm = 2rrdz + 2nzdr . (41) 


Because of equations (40) and (41) the change in energy dE due 
to a change in length dz is given by 


dE = 270 (4s — Yms) 2 + 2ryem(rdz + 2dr). (42) 
The volume of the cell-process is equal to 
V =x(?7? —a@*)z; (43) 
hence, since V = const., 
V 
ar = — dz. (44) 
2urz? 


In view of equation (44), equation (42) becomes 


dE = 220 (yes — yms) d2 + 2aven(T — )dz. (45) 


2urz 


As before, dH’ < 0, consequently since dz > 0 for elongation, we must 
have as an inequality equivalent to dE < 0 


V 
2na cs — ms fae mh) =. Us 
(Yes — Yms) + Zry¥em ( rae, (46) 
In view of equation (43), inequality (46) is equivalent to 
a? 2a 
LL aig (Yus.— Yoa) » (47) 


cm 


_, rts 
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in Figure 3 the expression 7 +- a?/r (the heavy line) is plotted against 
ry. Since 


! a? 2 
Aes feet ae area (48) 


we have a minimum value of r + a?/r atr =a. At r =a we have 
an ordinate of 2a. The dotted line a’ represents the value of 
20 (¥ms — Yes) /Yon LOY (¥ms — Yes) /Yem > 1. The dotted line b’ repre- 
sents the value of 20(yms — yes) /Yem for (yme — Yes)/Yon < 1. For 
(yms — Yes) /Yen < 1, inequality (47) cannot hold and thus elongation 
cannot occur. 

However, if (yms — yes) /Yem > 1, elongation can occur provided 
a<r<yr,. If r—ais of molecular dimensions, the basic assumptions 
are invalid. If r > r,, contraction (dz < 0) will result rather than 
elongation. Once elongation begins, it continues until either the end 
of the substrate fiber is reached, or until r—a approaches molecular 
dimensions. 

From equation (45) and (43), we obtain 

dE 


a ee 220 (ys a Vine) =~ MY em ( 
dz 


92 3 07 
) : (49) 


the 


Combining equations (49), (43), and (44) with equation (15) and 
the equation for the cross-sectional area in the present case, 


A, =a(7r?— a’), (50) 
gives 
dr a c( ie =the ) 
Fone GLEN ARS ted C8 1 icone | ere ate . (51) 
dt 3n Y va) 67 r 
Putting 
a Yom 
A=—— Maas Ads B= F) (52) 
37 (y ? Gy 
equation (51) becomes 
2q. ; 
ee Bde, (53) 
rP+—rt+a? 
B 


Integrated, equation (53) gives 
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A A 
r——log (7? +=7rt+ a) + 
Berra 


(54) 
A? 
— — 2a’ 2r + — 
B? 
— tan = aut aig 
4 —— AP 
J rR 
where 
A aries ; 
ies arg 47 ri G7 ia 
(55) 
A? 
— 2a? ary + = 
B? 
tan-? 
| A? A? 
4G AS 
B? E* 


Equation (54) is inconvenient for purposes of testing the theory. In 
theory we may obtain t as a function of z by introducing equation 
(43) into equation (54). Inspection of equation (51) shows that 
dr/dt is not equal to zero for r = a unless yms — yes = $7em- There- 
fore, we cannot expect the theory to have any physical significance 
over the entire range r 2 a. Another consideration limiting the the- 
ory is the fact that as r—a approaches molecular dimensions, the ba- 
sic assumptions are no longer valid. 

An equation giving dz/dt as a function of z may be obtained in 
the same way as was equation (51) except that r is eliminated in- 
stead of z. We have 


dz 2na TY om Vv a?\/nz 
Fs eRe (Yns'— Pee) 2? — 2 She f. . (56) 
dt 3nV 3nV ne  Va?naz+V 


Equation (56) might prove to be more convenient than an equation 
t = t(z) for purposes of testing. 


Ill 
Considering a frictional drag force of f. dynes per sq cm of Sg 


acting on the cell-process, the total force T pulling on the process is 
given by 


dE 
T =— — — 2naf,z. (57) 
dz 


ON ih \ ha a li ills 
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In order that elongation occur, this force must be positive, in other 
words, that 
dH 
——— 2aaf,z>0. (58) 
dz 
Hence, for dz > 0 


dE 


Combining equations (45) and (43) with inequality (59), we have 
@ (Yas — Yes Qa? 2a 
ser ais) Ze (reo); — (60) 
Yom r nV oe a a’) Vem 
as an inequality equivalent to inequality (59). 


Denoting 


a? 2af 
p(r) =r+—+ - 


Spee an eee (61) 
fe WV Yom (7? ao, a?) 


The general shape of the function y(7) is represented graphically in 
Figure 4. 
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FIGURE 4 


From equation (61) we obtain 


dy (r) went oa me Aafor ; (62) 
adr v2 wV Ym (7? a a?) 2 
Setting dy(r) /dr = 0, we have 
7? — 2 ez Aaf or (63) 
oe IV Voi To A 
Equation (63) becomes 
Aaf, 
(77 — a)? =r? —_—. (64) 
WV Vis 
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Denoting 
ee 
$= + ¥ Soe 2 (65) 


we obtain from equation (65) 
r?—a=sr. (66) 


Solving equation (66) gives 
So + V3? + 4a? 
——— = (67) 


Rejecting 8s, (—\/ 82 + 4a*) /2 since it is negative, we have as a positive 
value of r for which yw(7) has minimum 


8. + V8.2? + 4a? 
es (68) 


From Figure 4 and inequality (60) we see that if 
a ms 7 cs 
(y Ves) = 
Vem 
and if r lies between the two largest positive roots of the equation 


y(?), (69) 


eI = y(r) ; (70) 


then elongation is possible; the value of r decreases until it equals 
the smaller of the positive solutions of equation (70). This value will 
exceed 4. However, if r > a does not lie within this range, contrac- 
tion will occur. 
On the other hand, if 
a( CT cs) 
TS wlF), (71) 
Yom 
then inequality (60) cannot hold for any 7 > a and hence contraction 
rather than elongation will occur. 
Combining equations (15), (43), (45), and (58), we obtain 
dz 2na 227 om 


hoe me Yale 
at anv Yes) Z 3nV 


x a -f 7 Sareea deatne ® 2 _— 2n0fo2* : 
nz Ri V 3nV 
2nz4/ at + — 


NZ 


(72) 
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Integration of equation (72) would give z as a function of t, how- 
ever, the integration itself is a difficult task and for purposes of test- 
ing the theory, equation (72) might perhaps prove more convenient. 
than its integral. 

In all of these cases, we have failed to find an expression for 
dz/dt giving this derivative a constant value. The results of experi- 
mental work indicate that dz/dt is approximately constant at least 
through part of the total range of z. It is quite possible that other 
kinds of cel! substrate interfaces exist in cultures where these cells 
grow. A type of interface between cell and substrate in which con- 
tact between cell and substrate occurs only near the tip of the cell- 
process has been described. Furthermore, frictional forces may be of 
a different nature than those discussed here. The variability of the 
volume of the filopod may also play a role. Later theoretical investi- 
gations of this problem should be concerned with these equations. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Recent demonstration by the author has shown that the fundamental 
equations of the mathematical biophysics of the central nervous system 
can be considered as describing the behavior of very large numbers of 
neurons, of which each one follows discontinuous laws, such as discussed 
by W. S. McCulloch and W. Pitts. In that light some of the old problems 
are discussed. The comparative merits of the “microscopic” and “macro- 
scopic” approaches are discussed for the problem of the point to point 
correspondence between the retina and the cortex, with the number of 
connecting fibers much less than the number of cells. Some aspects of 
discrimination of intensities are also discussed. Finally, a few generali- 
zations of the McCulloch-Pitts treatment are suggested, and a nervous 
network is constructed which illustrates some aspects of the perception 
of numbers. 


In a previous paper (Rashevsky, 1945) we have shown that the 
usual equations of the mathematical biophysics of the central nervous 
system can be interpreted as a statistical macroscopic result of the 
interaction of a very large number of neurons, each following the 
laws now accepted as empirically established in neurophysiology. The 
“macroscopic” and the “microscopic” laws will from now on have to 
be used together, the former in those cases where a very large num- 
ber of neurons is involved, the latter for relatively small numbers. 
It is the purpose of this paper to consider briefly some problems which 
may involve the “microscopic” laws, for the description of which 
Boolean algebra is the adequate analytical tool (McCulloch and Pitts, 
1943; Householder and Landahl, 1945). 


I 

The problem of a point to point correspondence between the 
retina and some region of the occipital lobes has been treated from 
the “macroscopic” point of view by H. D. Landahl (1939). It does 
not seem to be possible on this basis to account for the fact that in 
certain parts the number of optical fibers is about one hundred times 
less than the number of cells which they connect. W. S. McCulloch 
(unpublished) has shown that by considering rather general temporal, 
patterns of discharge in a fiber, it is possible to construct a “micro- 
- scopic” mechanism which will give the desired effect. It is interest- 
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ing to consider as temporal patterns simple periodic discharges with 
different periods. 

Assuming with W. S. McCulloch and W. Pitts (1943) that all 
synaptic delays are approximately constant, being equal to about 0.5 
ms, and denoting this constant by 6, it is readily seen that a circuit 
as shown in Figure 1a will respond only when fiber I fires regularly 


eee a 


a b 
FIGURE 1 


with a frequency 1/6. Except where otherwise indicated, we sym- 
bolically assume (McCulloch and Pitts, 1943) that a simultaneous ex- 
citation of two terminal bulbs is necessary to make a neuron fire. In- 
hibitory synapses are symbolized by loops and we consider here only 
absolute inhibition. The circuit of Figure 1b will respond to a fre- 
quency 1/25. By interculating n internuncials, we obtain a circuit 
which responds to a frequency of 1/n6é. 

Hence, if a fiber discharges simultaneously with a large number 
of frequencies, all of the form 1/n6 , where n is a positive integer, and 
if this fiber branches off collaterals, each to a proper circuit, each cir- 
cuit will respond if, and only if, a corresponding frequency is pres- 
ent in the discharge of the fiber. 

Inasmuch as the usually observed steady frequencies in optic fi- 
bers are of the order of 20 per second or below, » must be no less 
than 100. If we wish the fiber to carry 100 frequencies of the form 
y = 1/n6 , then the highest one will be ymax = 1/1006 = 20 sec*, while 
the lowest will be »_,, = 1/2006 = 10 sec. There would be no diffi- 
culty in having one fiber per 100 retinal cells. A difficulty appears, 
however, in the necessity to assume that each retina cell connected 
to the fiber fires with a frequency 1/né, n being an integer, for if n 
is not an integer, no mechanism of the type shown in Figure 1 will 
respond to it. In fact, no circuit based on the assumption of a con- 
stant 6 will do that. 

Inasmuch as the retinal cells are seats of rather complex photo- 
chemical reactions, the possibility is open that the nature of the peri- 
pheral photosensitive processes provides for each cell a frequency of 
discharge of the form »y = 1/né. But at present this must be an ex- 
plicit additional assumption. 


II 
The mechanism for discrimination of intensities suggested by N. 
Rashevsky (1938b, 1940) and developed in detail by A. S. House- 
holder (1939, 1940a, b; Householder and Landahl, 1945) is based on 
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the macroscopic concept. Using, however, at the afferent end of a 
slowly adapting fiber such circuits as shown in Figure 1, it is seen 
that a stimulus of a given intensity S which results in a steady fre- 
quency »(S) may be relayed to different parts of the brain depending 
on its intensity. Empirically it is known that approximately 


y=logS. (1) 


But only frequencies of the form 1/né will produce a response at all. 
In order to have a mechanism ihat will account for perception of a 
continuous range of frequencies, we must consider the period of latent 
addition 7 , and thus obtain 


a : (2) 
Conta =). 
with - < 6. The smallest discriminable difference of frequencies is 
ps 1 1 1 1 
by =—_( = acs Je a . (3) 
2 7\ ni nmt1 +7n(n+1) 


For frequencies of the order of 200 sec"! or less, n > 10, and in- 
stead of equation (31), we have approximately 


tah 1 
dy = ——- —. (4) 
4) 2B ae 
From equations (1), (2), and (3) it follows that 
6S 
fae es) 10g" Ss) (5) 


which is contrary to experience. The old “macroscopic” theory 
(Householder and Landahl, 1945) is better in this case, although 
further modifications of the “microscopic” model may lead to differ- 
ent results. 

We may consider a group of rapidly adapting fibers so that the 
frequency in each fiber is independent of the intensity S of the stimu- 
lus. Let, however, due to differences in thresholds, the number of ex- 
cited fibers increase with S. The circuit shown on Figure 2 pro- 
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vides for an excitation of a different second order neuron for suffici- 
ently different values of S. But the relation between S and the J.N.D. 
6S is now determined solely by the distribution function of the thres- 
holds h;, which must be.chosen so as to give a correct relation be- 
tween 6S and:S. ; 

A distribution function of thresholds enters also into the “macro- 
scopic” theory of discrimination. All the work of A. S. Householder 
(1939, 1940a,.b; Householder and Landahl, 1945) indicates, however, 
that the gross, results of the theory do not depend: on the exact form 
of that distribution function-as they would in the present case. 


III 

We shall now discuss a mechanism for the discrimination of 
several temporally separated stimuli affecting the same end organ. 
On a “macroscopic” scale a somewhat similar problem has been dis- 
cussed by H. D. Landahl (1940). 

In the following, whenever logical symbolism is introduced, we 
shall, for typographical convenience, use the notations of D. Hilbert 
and W. Ackermann (1928) with the following two modifications: we 
shall use for equivalence the symbol ‘=” instead of “®”, and for 
conjunction a dot instead of “&”. 

The application of Boolean algebra to the theory of nervous nets 
is based on the circumstance that simple neuronic structures corre- 
spond to the functor S (McCulloch and Pitts, 1943), disjunction, con- 
junction, and negation, as illustrated by Figures la-1d of W. S. Mc- 
Culloch and W. Pitts, (1943). Just as any complex propositional 
function is built up by means of those four fundamental operators, 
so is any complex network built of the above.mentioned fundamental 
units. We shaJl introduce two new structures corresponding to the 
unlimited existential operator (Ex) N(x) meaning “there is an 2x 
for which N(«#) holds” and the limited existential operator (Hx) 
kN (a) meaning “there is an x amongst the integers from 0 to k for 
which N(x) holds.” That the structure corresponding to (Ex) is a’ 
closed circuit has been incidentally mentioned by W. S. McCulloch 
and W. Pitts (1943, p. 181, last line). Their illustration (1i) con- 
tains, however, a slip, in that the circle consists of a single neuron. 
This could not"be the case in the microscopic picture (Rashevsky, 
1945), where at least two neurons are necessary. 

Using 6 as the unit of time, we find the following: 

To the analytical expression 


N,(t) = (Ex) Ni (t—2x—1), (6) 


where « is an integer including zero, corresponds the circuit shown 
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in Figure 3a. To the expression 


N2(t) = (Ex) KN, (t—2—1), (7) 
Nii Ni2 Ni3 Nia 
Ny T 
No N, Ne 

a b 

N, nog peti No 
c d 
FIGURE 3 


corresponds the circuit of Figure 3b, with k neurons Nj; . 
A circuit, shown in Figure 3c, corresponds to 


Ni (t) = (Ex) Ni(t—2%—2), (8) 


and to the expression 
N.2(t) =N,i(t—1) - (Ex) N.(t—a4—1) tion) 


corresponds the circuit of Figure 3d. 

We now proceed to construct a network consisting of a single 
afferent N,, and of m efferents Ni, Nz --- N,,, and having the fol- 
lowing properties: Let N, fire at any irregular intervals, provided 
that all those intervals are integer multiples of 6, then when N, fires 
for the first time, N,, and only N,, responds; when it fires for the 
second time N., and only N., responds; when it fires for the j-th 
time N;, and only N; responds. Analytically this is expressed by 


N;(t) =N.(t—1) - (Ex) N,(t—a2—1) - (Ex) N.(t— 2-1): 
--» (Hz) N;.(t—«—1) (Ex) N;(t—a—1). 
But this is equivalent to the requirement Z 
N;(t) =N,(t—1) - (Ezy Nj.(t— 2 —1). (10) 
The corresponding circuits are shown in Figure 4. The neuron 


N, is assumed to have a threshold of 4, while the neurons of the rever- 
berating circuits and all other neurons Nj; have thresholds 2. 
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FIGURE 4 


If thus the same stimulus is applied to N, several times in suc- 
cession, each application elicits the response of a different N;. The 
activity of a given N; may release a set of motor reactions, verbal or 
otherwise, which will be different for each repetition of the same 
stimulus. We have here an elementary process of counting. 


If instead of the unlimited existential operator we use in expres- 
sion (10) the limited existential operator, a circuit is obtained simi- 
lar to that of Figure 4, but in which instead of circles, structures 
shown in Figure 3b are used. Such a circuit will have the property 
that repetitions of the same stimulus will be “counted” only when the 
intervals between the individual repetitions do not exceed a certain 
value determined by the value of k in (Ex)k. Speaking antropo- 
psychically, we may say that in this case we “forget’”’ the preceding 
repetitions if the intervals become too Jarge. This, perhaps, may 
correspond to actual situations. 


FIGURE. 5 
Consider now the circuit shown in Figure 5. Denoting by 
N(t, t+2) the statement that neuron N fires at the moments ¢ and 
t+2, the analytical expression for the circuit of Figure 5 is 


N.(t,+2) =N,(t—1)- 


a (11) 
(Ex) (N2(t—e—2) -Na(t—2))e 


In words: a single firing of N, results in a succession of two firings 
of N., provided N, has not fired in such a way before. We may now 
construct a circuit corresponding to 
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N;(t,¢+2) =N.(t¢—1) - (Ex) Nite 3, ¢ — = 1) 
- (Ez) No(t—x—38,t—2x—1) 
++: (Ex) Nja(t—2—3,t—a—1) So 
- (Ex) N;(t—ax—2,t—x—1). 


Returning again to one afferent N, and n efferents N;, we may 
conceive that N, is connected with each N; both by means of a 
circuit 3d, which is the fundamental unit used in constructing the 
circuit of Figure 4, and by means of circuits of Figure 5. In other 
words, we now construct a circuit which corresponds to both expres- 
sions (10) and (12). The circuit corresponding to expression (12) 
may, however, be inhibited by a permanently excited closed circuit 
C, which itself is inhibited by the single excitation of the n-th effer- 
ent N,. At the same time let the excitation of N, excite another 
closed circuit C, , which inhibits all circuits of the type of Figure 3a, 
at all the N;’s. Then up to 7 repetitions we shall have the situation 
described for Figure 4. But for the (n+1)st repetition of the stimu-~ 
lus, N, will fire twice, for the (n+2)nd repetition —N, will fire twice, 
and so forth. Putting 


(£2) Ni(t—2x) = [(E2) Ni(t—«x)] - 
. — oa” (18) 
-[(Ey) Ni(t—y)] -[e# #9], 


we have analytically 


N;(t) =N.(t—1) S&S (Ex) Ni(t—2—1)] 
-[(Ex)N;(¢—2—1)]- [7 <n]; 


i=j-1 


N,(t,t +2) =No(t—1) | > (Ex) N(t—2—8,t—2—1)| 


spb ax ot (14) 
*{ (£2) N;(¢—2—3,t—2—1)] 


SA tv ote oA YL, 


The circuit becomes too complicated to permit a practically useful 
drawing. ; 

The great disadvantage of the above mechanism is the require- 
ment that all stimuli ai N, be repeated at intervals which are integer 
multiples of 6, at least with an error not exceeding the period of 
latent addition. If this is not the case, the mechanism breaks down. 

The situation may be remedied in two ways. First, we may in- 
troduce some mechanism between N, and the rest of the system which 
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will provide for such a regulation of the intervals, even for stimuli 
spaced differently in time. Second, we may construct a circuit similar 
to the one represented in Figure 4, in which, however, pathways are 
substituted for fibers. In other words, we may consider a similar mech- 
anism on a “macroscopic” scale. The circuit then is considerably sim- 
plified since, due to the continuous excitation of the pathways, each 
neuron group N; needs to connect with only one pathway from each 
preceding Ni. (Figure 6). The group of neurons N, must be of a 


FIGURE 6 


mixed excitatory type (Rashevsky, 1938, 1940; Householder and Lan- ; 
dahl, 1945) characterized by 


A B 
A> Be. GoD Ste oe (15) 
a 6b 
For otherwise, a continuous stimulus to N, will produce a succes- 
Sive series of responses N,, N,2, N;, etc. Moreover, the duration of 
excitation of N, for a continuous stimulus must be less than the delays 
at the connections s; and s;' (Figure 6). Also, the saturation value 
E, of N. must be reached for very small values of the intensity S of 
the stimulus applied to N, so that the intensities E;’ of excitation of 
all the N;’ are independent of the intensity of the stimulus at N,. If, 
for simplicity, we assume that the constants of all the pathways con- 
verging on any connection s; are the same, and if we denote by E;” the 
intensity of excitation of each of those pathways, while #; denotes 
the intensity of excitation of the inhibiting pathways, then we must 
have 
AE; BE; AE; AE; BE; 
—— <h,<- + — — —.. (16) 
a b a a b 


With all the above requirements satisfied, the “macroscopic” net- 
work shown in Figure 6 will possess the necessary property. Various 
quantitative aspects may be studied by applying the fundamental 
equations of the mathematical biophysics of the central nervous sys- 
tem. 


eye ee et CAN 
= . xy 


N. RASHEVSKY 211 


It seems somewhat awkward to have to construct by means of 
Boolean algebra first a “microscopic circuit” and then obtain a sim- 
pler one by a transition to the “macroscopic” picture. We should ex- 
pect that a generalization of the application of Boolean algebra should 
be possible so as to permit its use for the construction of networks in 


which time relations are of a continuous, rather than of a quantized, 
nature. 
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A quantitative analysis of the immediate results of insulin and 
electric shock therapy given a scries of patients with schizophrenic re- 
actions leads to a definition of minimum standards of adequacy of 
treatment and suggests that the immediate outcome in adequately treat- 
ed patients depends on the duration of the illness. The experimental re- 
sults are rationalized by a formal theory which postulates a slow devel- 
opment of the illness, improvement to some extent after each treatment, 
and a slow relapse. 


Although a large number of papers on the results of shock ther- 
apy have appeared since 1936, little quantitative is known about it. 
It is to be expected a priovi that the results of treatments will depend 
on the number of shocks administered, as well as on the previous 
duration of the illness. The purpose of this paper is to make a step 
towards the study of the percentage P of recovery as a function of 
the number x of shocks and the duration 7; of the disease in months. 
Because of the small number of cases analysed, our results cannot 
claim anything more than a general indication of the type of existing 
relations. 

A preliminary analysis has been made of the results in 170 
courses of electric shock therapy and 185 courses of insulin shock 
therapy given in schizophrenic reactions at the Milwaukee County 
Hospital for Mental Diseases. The results were determined as of 
two to four weeks after the completion of therapy. 

Duration. The duration of the illness before treatment was de- 
termined from the history taken on admission. Where an estimate 
was given, as for example, six to eight months, the midpoint was 
taken as the time of onset. 

If a patient was discharged as recovered, then was readmitted 
and given shock therapy, the onset of illness was considered to be 
between the two admissions. If he had been discharged only as im- 
proved, he was assumed not to have recovered in the interval between 
admissions unless a very definite history to the contrary was ob- 
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tained; his illness was thus assumed to have begun before the first 
admission and to have continued until the time of treatment. 

Number of treatments. The number of treatments in the insulin 
group was defined as the number of days on which the patient be- 
came unconscious during the treatment; that is, only one coma could 
be experienced per insulin day. No distinction was made between 
long or short, deep or shallow comas. If no coma was produced, the 
patient was excluded from the series. The occurrence of a spontane- 
ous convulsion during the insulin day was not counted as an extra 
treatment. None of the patients received metrazol or electric shock 
therapy in conjunction with their insulin. 

The number of electric shock treatments was defined as the num- 
ber of times a current was passed through the patient, whether or 
not it produced a grand mal seizure. 

Insulin treatments ordinarily were given five times per week and 
electric shock treatments twice per week. Slight disturbances of this 
rhythm was disregarded in this study. A course was counted as ter- 
minated if more than thirteen days elapsed without treatment. All 
the cases discussed here received adequate treatment defined as having 
not less than thirty electric shock treatments, or -not less than fifty 
insulin comas. Values for nm greater than these limits are thus some- 
what arbitrary. 

Criteria of improvement. The degree of improvement was 
judged by the progress note written shortly after the termination of 
treatment. Recovery was taken to mean that the patient was free of 
symptoms and able to return to his former place in society. All lesser 
degrees of improvement were counted, in this study, as failures. 

Diagnosis. The diagnosis in all cases was schizophrenia. 

Sex. All the insulin patients were women. Both sexes were rep- 
resented in the electric shock group. Because in the large series of 
J. B. Ross and B. Malzberg (1939; Malzberg, 1943) there is shown no 
significant difference in the results of shock therapy in the two sexes, 
the electric shock results in this study were not shown according to 
sex, and the absence of men in the insulin group is not regarded as 
important. 

Results. In the total insulin group, 45% of the courses of treat- 
ment resulted in recovery; in the electric shock group the percentage 
of success was 41. The over-all results in the present series are of the 
same order of magnitude as those usually reported from other hos- 
pitals. It is therefore possible that the breakdown of the figures, 
which shows some interesting correlations, applies equally well to the 
reports from other institutions. For the RApduare groups the corre- 
sponding percentages are 76 and 59. 
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The results are shown graphically in Figures 1 and 2 for the 
cases which received adequate treatment. As said above, the small 
number of cases does not permit many conclusions. The general re- 
sult, however, seems to indicate that the function P(n, 7) increases 
with n, reaching an asymptotic value approximately at n = 30, and 
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that this asymptotic value is a decreasing function of 7. The indi- 
viduals having durations falling into arbitrary intervals were grouped, 
and to this group the mean duration 7 has been assigned. The inter- 
vals especially in the insulin data are large because the number of - 
cases is small. The number of cases in each group was as follows: 
for the insulin data from 7 = 2.8 to 7 = 63 the numbers are 37, 29, 
and 19, and for the electric shock data from short to long duration 
the numbers are 57, 34, 18, 32, and 29. 


Theoretical considerations. We shall make the following defini- 


tions: 
S =a measure of the state of the individual 
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H=S for average “normal” individual 

5 — decrement from average normal required to notice 
onset of abnormality 

A= decrement from average normal required to make 
diagnosis of return to normality (4 < 0) 

m = number of shock treatments, equally spaced in time 
and of equal intensity—first treatment origin of 
time scale 

i = value of m for which cure (S = H — A) takes piace 

+=duration of disease at initiation of therapy, time 
since S =H —o. 


Under normal circumstances one may assume that various events 
may act to change S,, but that homeostatic mechanisms influence the 
outcome. If Q(S, t), which may depend upon past events implicitly, 
is a measure of the effect of these events, S may be governed by a re- 
lation of the form 

dS 

dh ee (1) 
so that the normal value of S returns spontaneously toward H if 
the term Q does not exert too large an influence. If the effect of non- 
linearity of the homeostatic mechanism can be summarized by setting 
f(S) = (S—A)S, then equation (1) becomes 


ds 
dt 


=Q(S,t) —a(S—H) (S—A)S. (2) 


For Q=0,H,h, and 0 are equilibrium values of S, H and 0 being 
stable, h being unstable. Suppose now that S is in the neighborhood 
of H and that Q becomes large negatively and remains so for a long 
enough period of time. Then, according to equation (2), S will de- 
crease toward zero even after Q ceases to be negative, unless Q is too 
large positively. Consider only the case in which the average value 
of Q is negligible so that we may set Q = 0. Then roughly, for S <h, 
S decays exponentially according to A 


yes (H — 5) e-B(t+7) | (3) 


where b = aHh, and S = H — 6 for t= —-. This states that the 
condition of the individual spontaneously recedes from the normal, 
being, by definition, just noticeably abnormal (H—é) at t = — 1,7 
being the duration, and the time being measured from initiation of 
treatment. 


aa ee 
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The effect of the shock treatment is evidently to raise S eventual- 
ly. The data suggest that for the first few trials S may not be in- 
creased, but for the present we shall pass over this. In terms of equa- 
tion (2), this means that the treatment is equivalent to a positive Q 
over a period of time. Furthermore, S cannot continue to rise in- 
definitely so that the rate of rise of S with m must decrease for large 
n. This would occur if in addition to the terms on the right hand side 
of equation (2) there are terms which make for more rapid adjust- 
ment, but toward the value of S at some previous time, i.e., of the 
form — BS(t—7’). 

If ¢, is the constant time between treatments, then m = t/t, may 
be introduced as the variable. Furthermore, if the period of treat- 
ment is not too long so that S from equation (3) will not have been 
greatly altered, then for this period we will have 

dS 

—=akK—aS(t—7'), (4) 

dm 
where aK/t, is the effective Q due to treatment anda= ft,. If 7’ is 
larger than the period of treatment but too much sc, then S(t—7’) = 
S; from equation (3) with t set equal to zero, being approximately a 
constant. Hence from equations (8) and (4), ({< <7), 


S= (h—5d)e"+ K(1—e). (5) 
If m = n for which cure occurs, then 
H—<A=(H—d)e"+ K(1—e™”), (6) 


or if k= K/(H — 5), and e = (6 — A) /(H — 6) < 1, equation (6) 
becomes 
1—e7+e=—k(1—e™”). (7) 


Now k measures the relative effectiveness of the treatment to any 
individual. It is reasonable to expect that the logarithm of k is dis- 
tributed normally over the population. If p(z) is the normal curve 


1 Se 
for deviating z then we may write nf = log k/k ) , where o is the 
oT 


standard deviation and k& is the median value of the k’s in the popu- 
lation. From equation (7) we may then write for the per cent cured 
for those individuals with duration - and having had m treatments, 
the relation 
| P(n,7)= | p(z)dz, (8) 
1/o log k/k 

with k being given by equation (7). The curves in Figures 1 and 2 
illustrate the form of equation (8). 
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If, on the average, H — A < h, then more cases will relapse than 
remain cured since these will tend again toward the lower equilibrium. 
This will depend upon the value of Q which is generally unknown. 
Given the distribution of the h’s relative to H — A, and, disregarding 
the effect of Q, one can compute the percentage of relapses. At first 
glance this does not depend upon 7, but this is because we have not 
really taken into account the term 6 S(t — r+). Additional treatments 
after apparent cure would be beneficial according to this picture, es- 
pecially so for short durations. 

In Figures 1 and 2 are shown curves computed from equations 
(8) and (7) in which o = .41, b = .045, « = .67 for both sets of data, 
k = 1.57 and 1.69, a = .05 and .03 for the electric and insulin shock 
data respectively. If we take A = 0, and set H = 1, i.e., use arbitrary 
units with the normal] level unity, then 6 = .40, K = 1.94 and 1.01. 
The effectiveness per treatment, aK , is .047 and .030 in the two cases. 
From the point of view of the present theory o, 0, and « should be 
the same for both treatments (or any other) for equivalent groups 
since they are properties of the population. 

In the case of the insulin data the agreement with theory is not 
too unsatisfactory. It should be emphasized that the numbers in- 
volved are not large, while too large a range of values of + have had 
to be grouped together. In the case of the electric shock data, the 
agreement is less satisfactory. The principle deviation is for the 
short duration where the theory requires an appreciable difference 
between two and seven months’ duration, whereas none appears. How- 
ever, an appreciable difference could exist without showing up because 
the numbers involved are not large enough. 

The data for large + seem to reach an asymptotic value too soon 
to be accounted for by the theory. However, as pointed out above, the 
data for large n are particularly unreliable. But undoubtedly further 
modification of the theory will be required when additional data are 
available. ; 

Since the parameters are closely interconnected, the values given 
above are to be considered as very rough approximations. Further- 
more, the removal of simplifying approximations may change the 
shape of the curves somewhat while the values of the parameters may 
be changed considerably. 
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A NOTE ON THE MATHEMATICAL BIOPHYSICS OF 
CENTRAL EXCITATION AND INHIBITION 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


Some relations between the temporally macroscopic theory of cen- 
tral excitation and inhibition and the temporally microscopic theory of 
nervous nets are suggested. 


Consider a neuron ¢, which is acted upon by excitatory neurons, 
€;,°**+ , Ce, -*+, €, and inhibitory neurons Cp, ---, Cj, °°", Coq. Let 
6 be the least number of excitatory neurons, firing within a time 6, 
required to fire ¢,. The average frequency »v, with which c, fires will 
be given by [Householder and Landahl, 1945, chap. xv, equation (3) ] 


Yn — 6° (1—63>»;) ("v°), (1) 
where terms of degree greater than 6 have been neglected and where 
(*v°) is the sum of all possible terms of degree 4 in y of the p different 
y.s taken @ at a time. 

If p < 6, c, cannot fire. If p = 6, then », is proportional to the 
product of the »’s. If p= 6+ 1, there are p terms of degree 6. In 
this case, doubling any » practically doubles », . 

If the »’s are not of too different magnitude, one may write very 
roughly 


ee) 
n= 81-9 S95) es (2) 


p—6) 161’ 
where vy is some average value of the »’s, frequently between the 
arithmetic and geometric means. This can also be written approxi- 
mately as 

(6rv)(p—1)!Sr7 6°17? p! 


Vg rt 


bY jo 3 
(p — 0) 10! (p—8) 8!~” o 


peverally a better approximation than equation (2) if vy is the geo- 


metric mean and if p > 6+ 1. 
The quantity », can be considered as the stationary state fre- 


quency of c, and should thus be proportional to the quantity E of N. 
219 
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Rashevsky’s theory (Rashevsky, 1938). But then £ is linear with the 
E’s of the afferent fibers. Thus », should be linear in », (excitatory) 
and »; (inhibitory), whereas the first term on the right hand side of 
equation (3) contains »?". 

For the particular case in which only a subgroup m of the p 
fibers vary in magnitude, the others remaining constant in their ac- 
tivity, the effect of the coefficient »°* is reduced. In particular, if 
p—m=6—1,», is given by 

V_ = O°? G2 Veg Dy i MO Ty a. (4) 
v=1 
Thus for the case of several equivalent afferents acting upon a single 
neuron, approximate correspondence occurs for rather restricted con- 
ditions. 

We shall consider next some temporal relationships. If the fre- 
quency », determined from equation (2) depends upon time, due to 
variations in », and »;, so that »,. is the frequency at time sé, then 
the probability that ¢, fires for the first time in the interval 0 to ¢ is 
given by 


S=t/§ r=s-1 
P(t)= > I (1 = O¥ar) OV ne (5) 
s=1 Tr=1 
For constant »,, we have 
P(t) =1— (1— »,6)*, (6) 
For 6%, < <1, the above equation reduces to 
} P(t) =1—e™, (7) 


This suggests a possible identification of P(t) with the excitatory 
factor. If we set P(t) = 1/2 (or some arbitrary value), we can de- 
termine the value ¢, of t for which the probability is one-half that 
the first firing occurs prior to that time, so that t,, is the median re- 
action time given by 


(8) 


From the two-factor theory, the time to reaction will be given by 
the solution of a complex equation. If AF/a for each excitatory affer- 
ent is large compared with the threshold h, and a similar inequality 
holds for the inhibitory neurons, then one may write approximately 


h, 
SAE. — SBE, : (9) 


, 
m 
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for the median reaction time. 


But if », can be approximated by equation (4), the similarity be- 
tween equations (8) and (9) is evident since HE, « » and H; « »;. 
However, the latter expression involves a linear combination, the for- 
mer involves the sum. 
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A SUGGESTION FOR ANOTHER STATISTICAL INTERPRETA- 
TION OF THE FUNDAMENTAL EQUATIONS OF THE 
MATHEMATICAL BIOPHYSICS OF THE 
CENTRAL NERVOUS SYSTEM 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


In a preceding paper an interpretation of the « and j factors has 
been given in terms of an average effect of a large number of inter- 
neurons. In the present paper, a different interpretation is given in 
terms of the probability of a sufficient number of afferents to fire within 
the period of latent addition of the efferent. From this interpretation it 
follews that the old equations for « and j are only first linear approxi- 
mations to more complicated equations, the nature of which is suggested 
by this interpretation. 


In a previous paper (Rashevsky, 1945) we suggested a statisti- 
cal interpretation of the « and 7 factors. These were considered as 
measuring the number of self-circuited interneurons. The “synapse” 
of the mathematical biophysicist thus becomes a rather complicated 
structure, involving a very large number of interneurons. 

In the present note we outline a different approach to the prob- 
lem. 

Consider an efferent neuron whose threshold in terms of the 
minimum number of terminal bulbs which must be excited is n,. Let 
N afferent fibers synapse with this neuron, and consider for simplic- 
ity that they all have only one terminal bulb each on the neuron. Let 
n,<N. Let the average frequency vy of discharge be the same in all 
afferents. Let 6 be the period of latent addition. 

The probability for any one impulse to arrive at the neuron with- 
in a given interval 6 is equal to 6vy. The probability for n given fibers 
to bring impulses to the terminal bulbs within the interval is (dy)”. 
The probability P(n) for any n fibers to bring their impulses to the 
efferent neuron within the interval 6 is 


N! 


P (i) = (6r)” TITS 


(1) 


The smaller P(”), the longer we shall have to wait for such an 
event to happen. If the firing of the group of afferents begins at 
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t = 0, then the probable time ¢, of waiting is equal to 6 for P(n) = 1 
and to » for P(n) = 0. Hence, we may put approximately 
t) 


ar (n) 
Instead of N afferents, we may consider only one afferent, which 


connects with the efferent through N internuncials, so that the effer- 
ent is bombarded by irregular valleys of average frequency ». 


(2) 


ty 


Put 
dP (n) 
See a). (3) 
dn 
Equation (2) gives 
D) é 
P(n) =—; dP (n)=— — dt). (4) 
ty i 
Hence, from equations (3) and (4), 
6 dt, 
—=f(n). 5 
ee (5) 
But from equation (4) 
o2 
tS. 6) 
Seer ( 
hence 
P?(n) dt, 
aeris eee x 
pre f(n) (7) 


or, if we consider the number of simultaneously excited terminal 
bulbs as a function of the probable time t, , we have 


dn P*(n) 

dt, df (n) ~ 
The function P(n) and hence also P?(n) is finite for n = 0, and is 
zero for n = N + 1, being positive in that interval. The function 
f(m) is finite everywhere, provided (6v) is sufficiently small and N 


not too large, so that P(m) is monotonically decreasing with n . There- 
fore, in this case, the function 


(8) 


(9) 


is a positive monotonically decreasing function in the interval 
(0, N+1). 


As a first linear approximation we may put 


F(n) =A,—an. (10) 
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Since the value N + 1= A,/@, at which F (n) becomes zero, increases 
with N , we have approximately, with A, = A.N , 

F(n) =A.N—an, (11) 


or, since N + 1 « E, the total afferent excitation, (Rashevsky, 1945) 
therefore 


F(n) =AE—an, (12) 
and hence, from equation (8), 
ee Ty —an. (13) 
dt, 


Equation (13) is formally identical with the equations for « and 
j. Thus « is interpreted as the probable number of terminal bulbs ex- 
cited within the period of latent addition. The longer t,, the larger 
this probable number, but it never can exceed N. Hence if n, > N, 
then no matter how long we wait, the efferent neuron will never fire. 
This is the equivalent of the old requirement that h < AE/a, since 
nm, plays the role of h, and N « AE/a. 

Similarly, if we consider relative inhibition (McCulloch and Pitts, 
1943), 7 may be interpreted as the probable number of inhibitory syn- 
apses excited within the peried of latent addition. 

According to this interpretation, equation (13) is only a very 
rough approximation to the exact equation (8). The study of the 
latter thus leads us to a generalization of the present theory of cen- 
tral excitation and inhibition. 

The time t, at which x becomes equal to 7, is now also only an 
average value. The further development of this probabilistic theory 
will give us expressions for the fluctuations of t,, and hence of the 
reaction times. There are numerous data with which such expressions 
may be compared. 

The next important question now arises. Can we cient at least 
approximately, the relation E, = 6(e — 7 — h,) for the intensity of 
excitation # of the efferent fiber, at such times at which n > n,? With 
the special case, which we have so far considered, namely, N affer- 
ents synapsing with one efferent, this does not seem to be the case. 
For the quantity t, in equation (2) will also represent the average 
duration between two successsive events in which any 7 bulbs fire 
within the interval 6. For a time t > t, the most probable number 
n of bulbs firing within the interval 6 will be greater than 7, , but 
because of the all-or-none nature of the nervous discharge, this will 
not produce any changes in the latter. 

The picture of N afferents synapsing with only one efferent is, 
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however, rather artificial. It is much more plausible that on the effer- 
ent side we have to consider not a single fiber, but also a pathway 
consisting of N; fibers, all of different thresholds 7,(:) , and that each 
afferent fiber synapses with every efferent fiber. 

In this case, the efferent fiber with the lowest threshold , will, 
on the average, fire before all other fibers. This will happen when 
N= nN, , that is, at the time t,. For t > t,, n will exceed n,™ , the next 
lowest threshold, at a moment t,”, and so forth. After an infinite 
time, » will have reached its maximum value N, and all efferents 
whose thresholds are less than N will be firing, each with an average 
frequency »; = 1/t,. As EF, and hence N,, increases, the number of 
efferent fibers with thresholds 7,“ < N will also increase. Hence H#, 
will increase with 2, or with N. The increase will in genera! not be 
linear. The linear relation again is only a first approximation. 

If N is rather large, then the expression (1) for P(7) may have 
a maximum. In this case f(n) has a root at a point #, in the in- 
terval (0, N + 1), and therefore, F'(n) < 0 for n < %, becomes nega- 
tively infinite at that point, then changes sign and drops to zero for 
n—=N-+1. For the case of a constant N (viz. constant #) this would 
mean that the 7, t, curve (or ¢, t curve) will begin to rise only if it 
is already above n,, thus introducing a new sort of threshold. All of 
these possibilities indicate possible improvements in the present the- 
ory which are likely to cover a wider range of phenomena. 

That the present equations for « and 7 must be considered only 
as formal linear approximations to some more complicated equations 
based on physiological concepts has been pointed out previously (Ra- 
shevsky, 1940). The present paper is a step towards establishing 
these more exact and complicated equaticns. 

The author is indebted to Dr. H. D. Landahl for a discussion of 
the paper and for critical comments. 


This work has been aided in part by a grant from the Dr. Wallace 
C. and Clara A. Abbott Memorial Fund of the University of Chicago. 
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LETTER TO THE EDITOR 


THE HETERARCHY OF VALUES DETERMINED BY THE 
TOPOLOGY OF NERVOUS NETS 


Since you received the article on the above topic (Bull. Math. 
Biophysics, Vol. 7, No. 2), questions have arisen concerning the de- 
pendence of the anomaly on the topology of nets. It would be helpful 
if you could publish the following note: 

Given three dromes, each of which goes over one synaptic con- 
nection which is singly insufficient to fire its subsequent neuron but 
which may be reinforced from one other drome, an organism which 
fails to respond appetitively to any one of three sensory queues singly 
may respond to two by what appears to be a preference for one; and 
three such specicus choices may exhibit the circularity of the value 
anomaly. ; 

Consider three dromes—A, B, and C—so connected that no one 
sustains activity without summation from the afferent component of 
one other drome and let the net be such that B (and only B) neces- 
sarily contributes to A; similarly, C (and only C) to B, and A (and 
only A) to C. Presented with a stimulus a, b, or ¢ separately, there 
will be no response; but given any pair, a and b, or b and ¢, or ¢ and 
a, the organism will appropriate a, b, or c, respectively; and given 
a, b, and ¢, the organism will appropriate all three. Obviously the net 
resembles Figure 4 of the article in question except that the threshold 
of the afferent neurons is now such as to require impulses from the 
terminations of two axons, and that the heterodromic actions are sum- 
mative instead of inhibitory. The same topological considerations 
apply. The preference, whether or not it be a true choice, is deter- 
mined by a diadrome, which is no less a diadrome because its hetero- 
dromic connections are summative instead of inhibitory. 


WARREN S. MCCULLOCH, M.D. 
Illinois Neuropsychiatric Institute 
912 South Wood Street 
Chicago, Illinois 
June 20, 1945 


227 


INDEX TO VOLUME 7 


INDEX OF AUTHORS 


BRANSON, HERMAN. The Flow of a Viscous Fluid in an Elastic 
Tube: A Model of the Femoral Artery 


DANZIGER, LEWIS AND H. D. LANDAHL. Some Quantitative As- 
pects of Shock Therapy in Psychoses 


HILL, TERRELL L. AND LAURA E. HILL. Contribution to the The- 
ory of Discrimination Learning 


HOUSEHOLDER, ALSTON S. Muscular Dynamics and Muscular Ef- 
ficiency: I. The Isometric Length-Tension Diagram of 
Striated Skeletal Muscle 


HOUSEHOLDER, ALSTON S. Dynamics of Quadrupedal Locomotion 


KREUTZER, F. L. AND MANUEL F. MORALES. Some Nutritional and 
Excretional Interactions and the Growth of an Organ or 
Colony 


LANDAHL, H. D. Neural Mechanisms for the Concepts of Differ- 
ence and Similarity fF Perey 3 


LANDAHL, H. D. A Note on the Mathematical sae ee of Cen- 
tral Excitation and Inhibition : 5 


LANDAHL, H. D. AND LEWIS DANZIGER. Some Quantitative As- 
pects of Shock Therapy in Psychoses BEI Pe 


MCCULLOCH, WARREN S. A Heterarchy of Values Determined by 
the Topology of Nervous Nets Se eee a 


MORALES, MANUEL F. AND F.. L. KREUTZER. Some Nutritional and 
Excretional Interactions and the Growth of an Organ or 
Colony 


MORALES, MANUEL F. AND ROBERT E. SMITH. A Note on the 
Physiological Arrangement of Tissues ; mt eh 


MORALES, MANUEL F. AND ROBERT E. SMITH. The Physiological 
Factors Which Govern Inert Gas Exchange j 


OPATOWSKI, I. On the Form and Strength of Trees: Part III. 
The Secondary Branches and Their Relation to the Pri- 
mary Branches CER ep iewie ee oo hse Macs y ohne 

OPATOWSKI, I. Chain Processes and Their pee pee 
tions: Part I. General Theory ; : 


229 


PAGE 


181 


213 


107 


Or 


53 


15 


83 


219 


213 


47 


161 


230 INDEX OF AUTHORS 


RAsHEvsKy, N. A Problem in the Mathematical Biophysics of 
Blood Circulation: I re ta yi 
RASHEVSKY, N. A Problem in the Mathematical Biophysics of 
Blood Circulation II. Relation Between Pressure and 
Flow of a Viscous Fluid in an Elastic Distensible Tube 


RASHEVSKY, N. A Contribution to the Mathematical Biophysics 
of Visual Aesthetics a 5, Ae eee 


RASHEVSKY, N. Outline of a Mathematical Approach to the Can- 
cer Problem oe eee ee 


RASHEVSKY, N. The Mathematical Biophysics of Some Mental 
Phenomena é 


RASHEVSKY, N. Mathematical Biophysics of Abstraction and 
Logical Thinking ee HOOON Ee 

RASHEVSKY, N. A Reinterpretation of the Mathematical Bio- 
physics of the Central Nervous System in the ey of 
Neurophysiological Findings A a oa are 


RASHEVSKY, N. Some Remarks on the Boolean Algebra of Ner- 
vous Nets in Mathematical Biophysics Se ae 


RASHEVSKY, N. A Suggestion for Another Statistical Interpreta- 
tion of the Fundamental Equations of the Mathematical 
Biophysics of the Central Nervous System 

RUNGE, RICHARD. A Theory of Photosensitivity of Some Lower 
Animals ole ee ae cond eet ee 


RUNGE, RICHARD. A Contribution to the Mathematical Bio- 
physics of Cell Growth and Shapes: I : 


SMITH, ROBERT E. AND MANUEL F. MORALES. A Note on the 
Physiological Arrangement of Tissues 


SMITH, ROBERT E. AND MANUEL F. MORALES. The Physiological 
Factors Which Govern Inert Gas Exchange ee 


WILKINS, J. ERNEST, JR. The Differential Difference Equation 
for Epidemics , es hn ip id ee oo peas ee en 


PAGE 


25 


115 


133 


151 


203 


223 


59 


189 


47 


99 


149 


INDEX OF SUBJECTS 


231 


INDEX OF SUBJECTS 


Abstraction and logical thinking, 133 
Aesthetics, visual, 41 


Biophysical applications of chain 
processes, 161 

Blood circulation, 25, 35 

Boolean algebra of nervous nets, 203 


Cancer peomes mathematical approach 


69 
Cell cae and shapes, 189 
Central excitation and inhibition, 219 
Central nervous system, 151, 
Chain processes, biophysical applica- 
tions of, 161 
Concepts of difference and similarity, 
neural mechanisms for, 83 


Differential difference equation for 
epidemics, 149 

Discrimination learning, 107 

Dynamics of quadrupedal locomotion, 53 

Dynamics, muscular, and muscular ef- 
ficiency, 5 


Epidemics, differential difference equa- 
tion for, 149 


Flow of a viscous fluid in an elastic 
tube, 35, 181 
Form and strength of trees, 1 


Growth of cell, 189 
Growth of an organ or colony, 15 


Inert gas exchange, physiological fac- 
tors governing, 99 


Locomotion, quadrupedal, 53 
Logical thinking, 1383 


Mathematical approach to the cancer 
problem, 69 
Mental phenomena, 115 
Muscular dynamics and efficiency, 5 
Nervous nets 
boolean algebra of, 203 
topology of, 89, 227 : 
Neural mechanisms for concepts of dif- 
ferencé and similarity, 83 
Photosensitivity of some lower animals, 
59 
Physiological arrangement of: tissues, 
Physiological factors governing inert 
gas exchange, 99° 
Psychoses, shock therapy in, 2138 
Quadrupedal locomotion, 53 
Shock therapy in psychoses, 2138 
Tissues, physiological arrangement of, 
AT 


Topology of nervous nets, 89, 227 
Trees, form and strength of, 1 


Viscous fluid, flow of in an elastic tube, 


Visual aesthetics, Al 


ERRATA 


PAGE 


R. E. SMITH AND M. F. MORALES, On the ge of Ha Ele ae ete 


I. Fundamental Equations 


BOOK REVIEW 


ALSTON S. HOUSEHOLDER AND HERBERT D. LANDAHBL, Mathematical Beenie 
of the Central Nervous System (Reviewed by Herman Branson) 


SCOPE OF THE BULLETIN 

1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4, Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 


PREPARATION OF MANUSCRIPTS 


All manuscripts should be typewritten double space with margins 
14” on top, bottom and sides. Equations should be numbered on the 
right consecutively. Do not use such notations as “equation 2a” or, 
“equation 5” etc. References should be all given at the end of the 
paper, arranged in alphabetic order by authors and for each author 
chronologically, following exactly the style used in the Bulletin. In 
text, reference should be made by giving in parentheses the name of 
the author followed by the year of publication. In case of several pub- 
lications by the same author in the same year, use notations “1940a”, 
“1940b”, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
of the finished reproduction. Lettering should be made by means of a 
lettering guide. If desired by the authors, drawings may be made at 
the Editorial Office according to author’s instructions, at cost. 


